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ABSTRACT 


Two  computer  methods  for  the  design  of  radial  out-flow 
turning  vanes  are  described.  The  key  requirement  was  for 
vanes  with  simple  shapes  that  could  be  machined  easily  and 
inexpensively  in-house.  Two  design  approaches,  one  using 
circular  arcs  for  both  suction  and  pressure  sides  and  one  for 
wedge-arc  shaped  sections  were  developed. 

Requirements  for  the  general  design  procedure  arose  in  two 
exploratory  projects,  a  Rotary  Detonative  Wave  Engine  <RDWE) 
and  a  high  speed  Centrifugal  Diffuser  Test  Device  (CDTD> .  The 
RDWE  was  analyzed  first  and  found  to  be  potentially 
inefficient  in  its  original  configuration.  Preliminary  design 
for  a  high  speed  CDTD  flow  generator  was  completed, 
incorporating  wedge-arc  radial  outflow  vanes.  An  attempt  to 
analyze  the  flow  through  the  geometry  using  the  Eidelman- 
Godnov  Euler -(EGE^.  was  initiated. 
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NOMENCLATURE 


A»  throat,  area 

Ah  throat  area  for  High  Speed  CDTD 
Ai  throat  area  for  Low  Speed  CDTD 
a  throat  dimension 
b  width  teat  section 
D  out«ar  diameter 
e  energy  of  the  unit  volume 
L  length  of  blades 
m  mass  flow  rate 
P  power 

Pa  specific  power 
p  pressure 
R  radius  of  curvature 
ri  inner  radius 
r2  outer  radius 
s  blade  spacing 
T  torque 

U  velocity  component  corresponding  to  x  coordinate 
Un  velocity  component  normal  to  edge 
U2  wheel  rim  velocity 

V  velocity  component  corresponding  to  y  coordinate 
Vt  velocity  component  tangent  to  edge 
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V 9i  tangential  component  of  inlet  velocity 
V02  tangential  component  of  outlet  velocity 
W2  relative  exit  velocity 
Wx  work 

w'  rotational  velocity 
Z  number  of  blades 
a  efflux  angle 

B2  blade  exit  angle  with  respect  to  radial  direction 
82'  negative  blade  exit  angle  with  respect  to  radial 
direction 

r  specific  heat  ratio 
(  internal  energy 
*c  power  conversion  efficiency 
£  density 

Notations  used  in  Appendices  B  and  C  are  independently 


defined  therein 


IV.  FLOW  ANALYSIS 

A.  METHOD 

A  FORTRAN  computer  program  for  solving  numerically  the 
two-dimensional  Euler  equations  using  the  Godunov  method  was 
developed  at  the  Turbopropulaion  Laboratory  by  Dr.  S. 
Eidelman  [Ref.  12].  With  Dr.  Eidelman's  help  the  program 
was  to  be  applied  to  obtain  a  prediction  of  the  flowfield 
through  the  vane  passage.  An  abbreviated  description  will 
be  given  here  of  the  computational  method. 

The  unsteady  two-dimensional  Euler  equations  are  written 
as 

b  a  ♦  Bb  ♦  <b  c  =  0  C 7 ) 

b  t  bx  by 

where 

(O  PU  pU 

a  =  U,  b  =  p*pu2,c=*  PUV 

V  UV  p+  P  V2 

e  <e*p)U  (e+p) V 

and 

2  2 

e*  £  <  (  +  <U  ♦  V  >/2  >  *  energy  of  a  unit  of  volume 

(  * _ P  =  internal  energy 

<r  -Dp 

^  *  density 

U  and  V  are  the  velocity  components  corresponding  to  the  x 
and  y  coordinates 
p  =  pressure 
r  =>  specific  heat  ratio. 


25 


Modifications  were  aade  to  the  double  circular  arc 


design  to  achieve  convergency  at  the  exit.  The  first  step 
was  to  remove  the  pressure  line  and  use  the  camber  line  as 
the  pressure  surface.  The  assumption  was  made  that  the 
effects  associated  with  inlet  flow  incidence  resulting  from 
the  removal  of  the  pressure  line  would  be  negligible.  The 
results  indicated  that  the  passage  was  indeed  converging 
toward  the  exit.  The  second  step  in  an  attempt  to  assure 
tangency  and  convergency  at  the  exit  was  to  allow  for  a 
thickness  and  a  straight  section  at  the  leading  and  trailing 
edges.  The  straight  section  was  such  that  a  specified 
throat  dimension  was  achieved  at  the  trailing  edge  of  the 
adjacent  previous  vane.  The  resulting  program,  which  is 
presented  in  Appendix  C,  designs  a  wedge  shaped  arc  blade 
such  that  a  passage  which  converges  at  the  exit  is  realized 
with  a  specified  throat  dimension.  Using  the  method  a  high 
speed  CDTD  blading  design  was  obtained.  Figure  7  shows  the 
resulting  blade  design. 


literature  CRef.  8,9,10,11,  and  12  for  example!.  Yet  none 
was  found  which  dealt  with  the  required  case  of  radial- 
outflow. 

B.  DESIGN  METHOD 

The  expansion  process  in  a  radial  turbine  differs 
appreciably  from  that  in  an  axial  turbine  because  of  the 
radius  changes.  Vavra's  work  CRef.  51  involved  axial 
turbine  design.  The  work  presented  here  for  radial-outflow 
vane  design  followed  an  appreciation  of  Vavra's  methods  and 
possibly  retraced  steps  of  others  who  have  previously  sought 
the  simplest  geometries  for  machining  purposes. 

The  simplest  shape,  that  of  double  circular  arcs  was 
attempted  first.  Its  development  is  presented  in  Appendix 
B.  It  was  anticipated  that  with  the  requirements  of  radial 
inlet  flow  and  tangential  exit  outflow  the  double  circular 
arc  method  could  be  used  in  the  design  of  radial-outflow 
turning  vanes.  Once  the  method  was  developed  and  exercised 
it  was  clear  that  more  geometrical  constraints  were  required 
to  produce  a  nozzle  with  an  exit  which  was  both  nearly- 
tangential  and  convergent. 

The  double  circular  arc  program  generated  vanes  which 
appeared  to  converge  within  the  passage  and  then  diverge 
towards  the  exit  of  the  nozzle.  Therefore,  the  flow  would 
choke  within  the  nozzle  and  the  outlet  flow  conditions  could 
not  be  predicted  very  well. 


A.  BACKGROUND 


In  Ref.  5  Vavra  gave  a  procedure  to  eatabliah  a  blade 
profile  for  specified  valuea  of  throat  width  a.  spacing  s , 
trailing  edge  thickness  te»  and  inlet  flow  angle,  a.  Figure 
6  illustrates  these  features. 

Vavra' s  design  was  for  axial  turbine  blades  or  vanes.  A 
computer  prograa  was  available  to  construct  blading 
following  Vavra 's  Methods  [Ref.  5] .  The  method  used 
circular  arcs,  and  straight  lines.  Wedge  angles  and  leading 
and  trailing  edge  radii  were  input  parameters  to  determine 
a  flow  channel  with  suitable  acceleration.  No  comparable 
method  and  program  were  found  for  radial  outflow  vanes. 

In  1968  Goldman  and  Scullin  [Ref.  8]  designed  blades 
with  inlet  transition  arcs,  circular  arcs  and  outlet 
transition  arcs.  However,  these  blades  were  not  converging 
at  the  exit  and  were  not  designed  for  radial-outflow. 

Later  Goldman  and  Vanco  CRef.  9]  developed  a  computer 
program  for  the  design  of  sharped-edged  throat  supersonic 
nozzles  for  isentropic  flow.  Their  prograa  was  constructed 
for  the  radial-inflow  case.  The  Goldman  and  Vanco  design 
also  involved  the  requirement  of  contour  machining. 

There  is  a  substantial  amount  of  information  on  both 
radial-inflow  and  axial  turbine  blade  design  in  the 


<A*h/A«i>*  -  (6.27/38.3) a  *  .40.  (3) 

Using  this  scale  factor,  the  diaensions  calculated  for  the 


high  speed  sodel  are  shown  in  the  following  sketch. 


Figure  II. 2  Cross  Sectional  View  of  the  Test  Section 

From  equations  (2), 

a  z  *  A*  /  L  (6) 

therefore  for  a  diameter  of  8  inches  we  can  predict  a 
required  throat  dimension,  and  blade  spacing  using  equation 
<3>  for  a  selected  number  of  blades.  Results  are  given  in 
table  1.  The  requirement  then  is  to  design  suitable  turning 
vanes  to  produce  the  required  throat  area  and  spacing. 


Figure  II. 1  Description  of  Conditions  st  Throst 


The  blade  spacing,  s,  is  given  by 

a  ■  n  D  /  Z  (3) 

where  D  is  the  diaseter  at  the  vane  exit.  Froa  equations 
<2>  and  (3), 


A*  *  L  (a/s)  IT  D.  (4) 

In  Ref.  7,  the  following  conditions  were  calculated  for 

a  "full  scale"  high  speed  CDTD  design  to  have  a  test  section 

width  of  1  inch  at  25  inches  diaseter: 

Mia1.5  at  Bi«55°,  at  Ri*12.5  inches 
Vl»1392  ft/sec  at  To*520  »R,  giving 
V0i«114O  ft/sec  and  Vr0*1?57  ft/sec  and  for 
a  swirl  vane  length  of  5  inches,  Vr0«1B2  ft/sec 
and  Bo*82.34°.  An  increase  of  swirl  vane  length 
to  5.713  inches  gave  Vro>160  ft/sec  and  Bo*63.3° 


Using  the  above  data  for  the  "full  scale"  high  speed 
CDTD  in  equation  (4),  A*  ■  38.5  inches?.  Using  a  linear 
scaling  to  the  model  in  which  A*  *  6.27  inches?,  the  scsle 


factor  is 


C.  DESIGN  CALCULATIONS 

Vavra  [Ref.  53  showed  that  in  the  design  of  a  turbine 
cascade  the  efflux  angle  a  can  be  calculated  from  the 
geometry  of  the  cascade  and  that  it  is  independent  of  the 
incidence  angle  i  of  the  inlet  velocity,  aee  Figure  5. 

Vavra  also  noted  that  it  is  common  practice  in  steam  turbine 
design  to  use  the  relation 

a  *  arc  cos  <a/s)  (1) 

where  "a“  is  the  throat  opening  at  the  discharge  and  a  is 
the  blade  spacing.  Equation  <1>  applies  when  the  trailing 
edges  of  the  blades  are  of  negligible  thickness.  In 
practice  the  “wedge**  angle  at  the  trailing  edge  and  trailing 
edge  bluntneas  modifies  this  dependence. 

With  the  limitation  of  a  two  minute  runtime  and  the 
available  mass  of  stored  air.  Demo  [Ref.  63  found  that  for 
the  supersonic  cascade  rig  in  building  230  an  area  of  6.27 
square  inches  was  needed  at  the  nozzle  throat  where  the  flow 
would  choke.  This  same  area  was  taken  for  the  total  outlet 
throat  area  of  the  swirl  vanes  for  the  high  speed  CDTD. 

Referring  to  the  following  sketch,  the  area  at  the  exit 
where  the  flow  will  choke.  A*,  is  given  by, 

A»  •  L  a  Z  (2) 

where  L  is  the  vane  length,  a  is  the  throat  width  and  Z  ia 
the  number  of  vanes. 


required  across  the  swirl  generator  in  the  high  speed 
version,  turbine  like  turning  vanes  were  selected  rather 
than  the  sheet  metal  nozzle  construction  used  in  the  low 
speed  device. 

B.  APPARATUS  DESCRIPTION 

Figure  3  illustrates  the  overall  arrangement  of  the 
proposed  CDTD  design.  The  swirl  vanes  for  the  high  speed 
CDTD  are  sandwiched  between  two  circular  plexiglas  disks. 
Plexiglas  disks  were  used  to  allow  visualization  of  the  flow 
through  the  vanes.  The  disks  are  mounted  on  an  aluminum 
center  pipe  shown  in  Figure  3.  The  pipe  provides  double 
axial  in-flow.  Following  the  work  performed  by  Feieresien 
[Ref.  4]  a  transition  section  to  turn  the  air  from  the  axial 
to  the  radial  direction  was  added  in  the  form  of  a  plate  on 
a  tie-rod  on  the  pipe  centerline.  The  design  of  the 
transition  section  is  not  complex  if  the  area  of  the  radial 
outlet  is  smaller  than  that  of  the  total  axial  inlet.  An 
acceleration  through  the  transition  section  will  then  occur, 
leading  to  generally  favorable  pressure  gradients.  These 
consideration  required  a  double-entry  supply  arrangement  for 
the  CDTD  model . 

The  piping  is  attached  to  the  air  supply  system.  Figure 
4  shows  the  configuration  of  the  blades  on  the  plexiglas 
disks  and  in  relationship  to  the  supply  pipe. 


II.  APPARATUS  DESIGN 


Factors  considered  in  the  design  of  the  high  speed  CDTD 
fell  into  two  categories;  namely,  limitations  due  to  the 
available  air  supply  system  and  piping  and  the  requirement 
that  the  test  section  have  dimensions  appropriately  scaled 
for  the  diffusers  required  to  be  tested. 

A.  DESIGN  CONSIDERATIONS 

The  high  speed  CDTD  was  designed  to  operate  as  an 
intermittent  blow-down  wind  tunnel.  In  order  to  use  the  air 
supply  system  available  in  building  230  at  the  Naval 
Postgraduate  School,  the  test  section  flow  area  was  dictated 
by  the  available  mass  of  stored  air  and  the  minimum  desired 
runtime  of  two  minutes.  Using  proportions  similar  to  those 
used  in  the  low  speed  CDTD  <see  Figure  2>,  the  radius  of  the 
swirl  generator  became  four  inches  while  the  radius  at  the 
teat  section  was  five  inches.  The  method  used  for 
calculating  these  dimensions  is  given  in  the  following 
section.  The  difference  between  the  swirl  generator  outer 
radius  and  the  test  section  radius  allowed  for  mixing  of  the 
wakes  from  the  blading  (or  nozzles)  used  in  the  swirl 
generator.  The  swirl  generator  was  limited  to  a  length  of 
two  inches  to  meet  the  scaling  requirements  set  by  the  test 
section.  Because  of  the  larger  pressure  differential 


apparatus  are  described.  General  methods  and  computer 
programs  for  the  design  of  double  circular  arc  and  wedge  arc 
blading  for  radial-outflow  vanes  are  described  in  detail  in 
appendices.  The  report  also  discusses  a  method  of 
predicting  the  flowfield  through  a  selected  passage  using  an 
Euler  code  based  on  the  Godunov  method  of  solution. 


development  for  viscous,  three  dimensional,  transonic  flow 


(a  regime  containing  both  supersonic  and  subsonic  regions), 
with  adverse  static  pressure  gradients.  Sucn  methods  are 
required  to  design  optimum  radial  diffusers  for  transonic 
conditions.  Equally,  a  high  speed  version  of  the  CDTD  was 
required  in  order  to  examine  such  diffuser  flows  under 
controlled  conditions  and  to  verify  the  emerging  codes. 
Again,  nearly  tangential  swirling  outflow  was  required,  this 
time  in  steady  flow. 

The  main  interest  here  was  therefore  to  design  vane 
shapes  and  the  apparatus  for  a  high  speed  CDTD.  The  design 
would  closely  parallel  the  general  arrangement  of  the  low 
speed  CDTD.  The  high  speed  CDTD,  had  the  unusual 
requirement  of  radial-outflow  with  convergency  at  the  rim 
exit.  The  limitation  on  the  design  was  primarily  that  the 
hardware  could  be  machined  easily  and  inexpensively  in- 
house.  Computer  programs  for  the  design  of  supersonic 
nozzles  for  radial-inflow  turbomachinea  were  found  to  exist, 
but  most  were  highly  involved  and  would  need  extensive 
machining  processes  to  generate  the  blade  shapes. 

Particular  designs  are  discussed  herein.  A  literature 
search  indicated  that  computer  programs  for  radial-outflow 
devices  using  simple  vanes  for  two-dimensional  isentropic 
flow  of  a  perfect  gas  were  not  available. 

In  the  present  report,  the  design  of  the  blading  for  a 
high  speed  CDTD  model  and  the  general  arrangement  of  the 
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The  question  then  was,  what  would  be  the  best  design 
for  the  rotor  passage  to  maximize  the  torque  ispulse?  In  an 
analysis  included  in  Appendix  A  it  was  shown  that  tangential 
flow  at  the  ris  of  the  RDWE  would  yield  the  highest  specific 
power  and  highest  power  conversion  efficiency.  A  design  was 
required  therefore  of  nozzles  which  would  direct  the  outflow 
alsost  tangentially.  The  nozzles  would  operate  in  an 
unsteady  sode. 

The  second  project  which  Motivated  the  present  study 
involved  the  developaent  of  a  test  rig  for  radial 
diffusers.  Such  diffusers  are  used  in  centrifugal 
cospressors  to  convert  the  kinetic  energy  of  the  flow  from 
the  rotor  into  static  pressure.  Erwin  [Ref.  21  proposed  and 
designed  a  so-called  low  speed  Centrifugal  Diffuser  Test 
Device,  (CDTD) .  The  CDTD  simulated  the  flow  of  air 
delivered  to  the  diffuser  from  the  compressor  rotor  using 
static  swirling  vanes.  The  facility  would  allow  the 
evaluation  of  proposed  new  diffuser  geometries  and  permit 
validation  of  computer  analyses  for  diffusers  operating  in 
the  fully  subsonic  flow  regime.  Figure  2  shows  a  schematic 
of  the  low  speed  CDTD.  Vidoe  [Ref.  33  presented  results 
from  the  flow  generator  used  in  the  CDTD  and  performed  an 
analysis  of  the  internal  flow.  As  mentioned  by  Vidos,  the 
design  of  centrifugal  diffusers  is  presently  based  on 
experimental  results  for  two  dimensional  and  conical 
diffusers.  Numerical  methods  are  currently  under 


I .  INTRODUCTION 


The  Motivation  for  the  work  reported  here  caae  fro*  two 
exploratory  projects,  a  Rotary  Detonative  Wave  Engine, 
<RDWE>,  described  in  Ref.  1  and  a  high  speed  Centrifugal 
Diffuser  Test  Device  (CDTD),  described  in  Ref.  2.  Both 
projects  required  the  geoaetrical  design  of  an  axisyMjnetric 
generator  for  nearly-tangential  transonic  outflow. 

To  iaprove  gas  turbine  engine  perforaance  an  increase 
in  cycle  pressure  ratio  is  needed.  One  alternative  to 
Mechanical  coapression  is  to  supply  the  coapression 
thermally  by  detonative  coabustion.  If  the  pressure  rise 
produced  in  this  Manner  is  sufficient  for  the  engine,  no 
coapressor  and  associated  drive  turbine  are  required  and  a 
siapler  aore  efficient  engine  Might  result. 

As  described  by  Honks  CRef.  13  ,  a  proposed  geometry 
for  such  an  engine  aight  include  a  rotor  with  outflow 
passages  which  are  turned  at  the  exit  to  produce  nearly 
tangential  flow.  Figure  1  illustrates  the  proposed  RDWE. 
Honks  atteapted  to  measure  the  torque  on  a  static  simulated 
rotor  as  a  detonation  wave  propagated  froa  the  center 
outward.  A  reaction  torque  was  produced  from  the  nearly- 
tangential  expulsion  of  gases  froa  the  ria. 
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Equation  (7)  is  nondimensional ized  using  arbitrary 
reference  quantities,  and  the  Godunov  method  [Ref.  121  is 
applied  to  compute  at  each  time  interval  the  fluxes  at  the 
boundaries  of  the  cells  of  a  grid  established  for  the 
problem  under  consideration.  The  steady  state  solution  is 
found  when  no  further  change  occurs  between  successive  time 
steps. 

The  computational  domain  is  shown  for  a  typical  cascade 
problem  in  Figure  8  and  the  grid  and  cell  notation  are 
defined  in  Figure  9.  Referring  to  Figure  8,  for  the 
tangency  condition  at  the  wall  of  the  passage  a  solid  wall 
boundary  condition  is  imposed  over  the  segments  3-5  and  4- 
5.  Boundary  conditions  for  segments  1-3,  5-7,  2-4  and  6-8 
are  for  periodicity  of  the  flow  with  respect  to  the  y  axis. 

The  parameters  held  constant  at  the  inlet  are  the  x 
components  of  velocity,  U,  the  flow  angle,  9,  the  gas 
stagnation  enthalpy,  H,  and  the  entropy,  S.  The  velocity, 
density  and  pressure  can  be  calculated  from  these  quantities 
easily  as  Vin,  Pin»  and  pin*  At  the  outlet  boundary  only 
pressure  is  defined  and  other  flow  parameters  are  calculated 
in  the  solution.  In  the  solution  procedure,  Uout,  Vout,  and 
P  out  are  set  equal  to  the  values  of  U,  V  and  p  one  point 
ahead  of  the  outlet  boundary. 

A  subroutine  (program)  constructs  a  non-orthogonal  grid 
covering  the  computational  domain  as  shown  in  Figure  9.  The 
flow  parameters  related  to  the  center  of  the  cell  have  a 
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fractional  index:  k+&.  The  cell  boundaries  have  one 

fractional  and  one  integer  index:  3,  k*  H  or  m+!4,n.  The 
paraaeters  at  the  tine  t  will  have  subscript  indices  and  at 
tine  t*  St  the  parameters  will  have  superscript  indices. 

The  Euler  equation  is  approximated  to  first  order  by  the 
equation  in  table  2.  Solution  to  the  itiemann  problem  is 
initiated  assuning  that  the  left  state  is  that  located  at 
j-l^k+h  and  the  right  state  is  that  at  )4»k+ >4. 

The  solution  to  the  Riemann  problem  will  give 

<R,P,Un,vt,>3fk+Ji 

Un  3  the  velocity  component  normal  to 
the  edge  <j,k*l>  -<j,k) 

Vt  =  the  velocity  component  tangent  to 
the  edge  <j,k+l>  -<j,k> 

The  capital  letters  R  <rho>,  U,  V,  P  are  the  parameters 
calculated  at  the  edges  of  the  zones. 

The  component  of  the  velocity  is  then  transformed  back 
to  the  cartesian  coordinate  system  and  a  final  value  of  the 
parameter  on  the  cell  edge  is  received 

(R,P,Un»B)  ^ 


B.  RESULTS 

The  program  (EGE)  was  run  on  the  IBM  370-3033  at  the 
Naval  Postgraduate  School.  Figure  10  illustrates  the  grid 
covering  the  passage  for  the  wedge-arc  blading  shown  in 
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Figure  7.  An  enlargement  of  the  exit  area  is  shown  in 
Figure  11. 

In  order  to  apply  the  code  to  the  radial  geometry,  some 
modifications  were  necessary.  A  transformation  of 
coordinates  was  made  to  conveniently  accommodate  the  radial 
flow  inlet  boundary  condition.  Since  the  program  remained 
basically  two-dimensional,  the  outlet  flow  was  treated  as 
being  periodic  over  the  blade  space  at  the  exit  which  then 
did  not  change  with  radius.  Thus  the  solution  found  will  be 
valid  for  the  flow  within  the  rim  and  approximately  true 
just  outside  the  rim.  For  this  reason  the  computational 
domain  in  Figure  11  was  not  extended  radially  outwards 
beyond  one  blade  space. 

Results  could  not  be  obtained  for  the  flow  through  the 
passage  within  the  time  available.  Some  dif f iculities  were 
encountered  by  the  program  in  the  region  of  the  throat. 

Some  modifications  of  the  grid  must  be  made  before  a 
converged,  steady  flow  solution  can  be  obtained. 
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V.  SUMMARY  AND  DISCUSSION 


Two  exploratory  projects  at  the  Turbopropulaion 
Laboratory,  1).  to  test  the  concept  of  detonatlve  combustion 
as  an  alternative  to  mechanical  compression  in  gas  turbines 
and  2) .  to  design  a  teat  rig  for  radial  diffusers  used  in 
centrifugal  compressors,  were  examined.  The  projects  were 
similar  in  nature  since  they  involved  the  generation  of 
nearly-tangential  radial  outflow  and  each  required  a  design  of 
turning  vanes  for  a  teat  apparatus.  The  overriding  limitation 
for  either  project  was  that  any  apparatus  developed  would  have 
to  be  easily  and  inexpensively  manufactured  in-house. 

In  the  Rotary  Detonative  Wave  Engine,  "RDWE" ,  project, 
turning  vanes  were  required  such  that  upon  detonation  the 
maximum  possible  torque  impulse  would  be  obtained.  The  best 
design  possible  to  yield  maximum  power  hence  maximum  torque 
was  to  produce  turning  vanes  that  exit  nearly-tangential. 

Using  Euler's  Turbine  equation  a  relationship  was  developed 
relating  exit  angle  to  specific  power  and  efficiency. 

The  results  for  the  unsteady  expansion  process  showed 
that  the  value  of  the  relative  flow  angle  (82*  established  the 
limiting  conversion  efficiency  and  specific  power.  If  the 
relative  velocity  increases  beyond  the  limiting  point  the 
specific  power  generated  will  be  absorbed  rather  than 
generated  and  efficiency  will  decrease.  High  efficiency 


occurs  in  large  part,  only  when  the  specific  power  is  low  or  at 
the  Halting  point.  In  general  the  conversion  efficiency  will 
vary  during  the  expansion  process,  but  will  always  be  poor. 

This  result  implies  that  the  concept  of  the  RDWE  can  not 
be  efficient  aa  a  atand  alone  rotor.  Other  methods  would  have 
to  be  employed  to  increase  efficiency,  such  aa  the  utilization 
of  additional  components,  particularly  those  which  would  allow 
nearly  constant  flow  by  sequencing  the  detonation  wave  with 
respect  to  partial  porta. 

The  second  project  involved  the  development  of  a  test  rig 
for  a  high  speed  Centrifugal  Diffuser  Test  Device,  "CDTD". 

The  project  would  cloaely  parallel  an  in-house  concept  of  a 
low  speed  CDTD.  The  low  speed  version  simulates  the  flow  of 
air  delivered  to  a  diffuser  aa  a  steady  process.  The  high 
speed  version  simulates  in  a  steady  process,  the  transonic 
Mach  numbers  characteristic  of  advanced  centrifugal 
compressors . 

For  the  high  speed  CDTD  a  passage  would  be  required  which 
would  cause  the  flow  to  converge  and  choke  at  the  exit  of  the 
passage.  Because  the  concept  of  the  RDWE  was  questionable  the 
main  interest  tuned  to  designing  vane  shapes  and  the  apparatus 
for  a  high  speed  CDTD. 

Two  methods  and  their  respective  computer  programs  were 
developed.  The  first  method  involved  the  development  of 
double  circular  arc  vane  shapes.  Given  the  outer  radius,  the 
inner  radius  and  the  number  of  blades  the  radii  of  the  suction 


and  pressure  sides  were  determined.  The  design  emphasis  was 
on  acheiving  a  tangent  suction  surface  at  the  exit  and  a 
radial  camber  line  at  the  inner  radius.  The  concurrent 
condition  of  a  tangent  camber  line  at  the  exit  gave  the 
solution  of  a  blade  of  zero  thickness. 

The  geometrical  description  resulting  in  a  tangent  camber 
line  at  the  exit,  resulted  in  the  sine  of  the  camber  angle  0 
being  equal  to  the  cosine  of  the  subtended  blade  angle  y>  . 
Following  manipulations  to  relate  the  value  of  to  just  the 
radius  ratio  (the  only  specified  parameter),  an  equation  of 
fourth  order  in  sin  f*  resulted.  Using  Newton's  method  of 
interation  the  first  (smaller)  root  was  found.  Using  this 
value  of  'I'  a  blade  shape  was  generated.  It  was  clear  that 
the  auction  side  was  tangent  but  the  efflux  angle  a  was  not 
zero.  An  efflux  angle  of  zero  would  indicate  that  the  camber 
line  was  also  tangent  to  rim.  A  graph  of  the  fourth  order 
function  indicated  that  two  roots  were  possible  for  each 
radius  ratio.  Only  the  first  root  was  solved  for.  It  is 
probable  that  the  second  root  would  yield  the  tangency 
condition  of  both  camber  line  and  suction  aide.  The  selected 
root  was  shown  however  to  give  sin  0  *  cos  “f"  . 

It  was  realized  that  convergence  at  the  exit  could  not  be 
guaranteed  with  the  double  circular  arc  blade  design. 
Modifications  were  applied  to  the  design  to  cause  the  passage 
to  converge.  Improved  convergence  resulted  when  the  original 
camber  line  and  the  suction  line  were  used  to  describe  the 


blade.  The  aiaple  double  circular  arc  deeign  required  curved 
or  thin  trailing  edges.  With  the  addition  of  a  finite 
thickness  at  the  leading  and  trailing  edge,  a  straight  line 
segment,  and  a  wedge  angle  all  of  which  would  better  direct 
the  flow  in  the  tangent  direction,  practical  blade  shapes  with 
convergent  passage  were  acheived. 

Geometrical  equations  were  developed  to  yield  the  proper 
blade  design.  The  exit  throat  dimension  was  fixed  in  the 
design  based  on  the  limitations  of  the  available  air  supply 
and  baaed  on  dimensions  scaled  from  a  high  speed  CDTD  design. 
It  was  anticipated  that  the  flow  would  choke  at  the  specified 
throat  dimension  between  adjacent  auction  and  pressure  sides 
and  the  straight  line  segment  would  intersect  the  exit  circle 
which  was  tangent  to  the  auction  side.  This  did  not  occur. 

The  length  of  the  straight  line  segment  varied  between  1  and  5 
percent  of  the  value  needed  for  this  condition.  It  was  not 
required  that  the  exit  circle  be  tangent  to  the  pressure  side. 

From  the  geometry  and  the  specified  parameters  a  blade 
was  developed  giving  convergency  at  the  exit.  The  camber 
angle  0  and  the  subtended  blade  angle  ^  were  specified  from 
experience  obtained  in  the  development  of  the  double  circular 
arc  design.  The  chosen  angles  0  and  ^  gave  an  efflux  angle  a 
between  8  and  9  degrees.  To  achieve  the  perfect  condition  cf 
tangential  flow  angle  oc  would  be  required  to  be  equal  to 
zero.  Based  on  current  technology  and  manufacturing  abilities 
an  angle  in  the  area  of  8  or  9  degrees  would  give  as  close  to 
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tangential  flow  as  is  possible.  It  is  noted  that  the  angle  9 
is  related  uniquely  to  the  angle  J'  .  Once  the  radius  ratio  is 
given . 

The  coordinates  of  the  designed  passage  were  then  input 
into  the  FORTRAN  program  (EGE>  to  predict  analytically  the 
flowfield  through  the  vane  passage,  by  solving  the  Euler 
equations.  The  program  (EGE)  was  designed  to  run  on  the  IBM 
370-3033.  The  program  successfully  constructed  a  non- 
orthogonal  grid  covering  the  domain,  however  results  could  not 
be  obtained  for  the  flowfield  through  the  passage  in  the 
initial  attempts.  It 

is  thought  that  the  difficulty  lies  in  the  curvature  and 
density  of  the  grid  in  the  region  of  the  exit  throat.  Since 
the  EGE  program  has  not  previously  been  applied  to  a  radial 
outflow  problem  with  severe  curvatures  the  task  of  analyzing 
the  present  flow  with  the  EGE  code  is  now  thought  to  require 
the  concentrated  effort  of  the  code's  author. 

In  contrast,  the  CDTD  vane  design  can  be  manufactured 
without  difficulty  and  the  flow  from  the  turning  vanes  can  be 
measured  experimentally.  It  is  suggested  that  such  an 
experiment  would  provide  a  severe  test  of  the  ability  of  the 
EGE  code  to  simulate  transonic  radial  outflow  containing 


embedded  waves 
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VI .  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  following  concluaiona  are  offered: 

The  Rotary  Detonation  Wave  Energy  concept  in  ita  original 
configuration  ia  potentially  inefficient.  An  increase  in 
efficiency  would  require  the  addition  of  a  component  which 
could  recover  the  unateady  kinetic  energy  of  the  flow  leaving 
the  rotor . 

A  program  for  designing  radial  outflow  vanes  using  double 
circular  arc  blades  waa  demonstrated  successfully.  While  the 
method  may  be  useful  in  other  contexts,  it  waa  found  to  be 
unable  to  produce  converging  passage  shapes  with  a  near- 
tangential  outlet  flow  direction. 

A  program  for  designing  radial  outflow  vanes  using  wedge- 
arc  blade  shapes  waa  uaed  successfully  to  produce  converging 
passage  shapes  with  prescribed  throat  widths  and  with  a  near- 
tangential  outlet  flow  direction. 

Both  blade  shapes  can  be  produced  relatively  inexpensively 
from  standard  metal  atock. 

A  preliminary  design  of  a  test  model  of  a  high  speed 
Centrifugal  Diffuser  Test  Device  was  produced. 
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B 


RECOMMENDATIONS 


The  following  recommendations  and  areas  of  further 
investigations  are  suggested: 

An  arrangement  of  partial  porting  and  following  power 
turbine  for  the  Rotary  Detonation  Wave  Engine  should  be 
examined. 

The  program  for  double  circular  arc  blades  should  be 
modiied  to  seek  both  roots  of  the  condition  for  camber  line 
tangency,  and  to  offer  the  selection  of  any  value  of  ^  between 
them.  This  will  result  in  the  ability  to  examine  any  value  of 
yj  within  the  range  of  any  practical  interest. 

In  the  case  of  wedge-arc  blades,  the  possibility  of 
solving  for  the  particular  case  in  which  the  circle  at  the 
passage  outlet  is  also  tangent  to  the  pressure  side  should  be 
examined.  A  procedure  should  also  be  introduced  to  calculate 
the  progression  of  the  passage  area  taking  into  account  the 
increasing  radius. 

The  proposed  blading  geometry  should  be  constructed  and 
the  outlet  flowfield  examined  experimently . 

Work  should  be  continued  to  obtain  a  computational 
prediction  of  the  flow  through  the  passage.  The  model  test 
could  then  serve  as  a  verification  of  the  code. 
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TABLE  1.  Required  Design  Dimensions  for  Convergent. 
Throat  in  High  Speed  CDTD. 


NUMBER  OF  BLADES 
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TABLE  2.  First  Order  Solution  of 
Euler  Equations 
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hV2  =  h  U22  c  -2  W-?  sin  O2 '  +  1  +  _W2.  2  3  A<13> 

U2  U2 

Substituting  equation  A(10>  and  A(13)  into  A(9> 

"c  =  P/» 

P/a  ♦  )jV22 

=  U^2  <W2/U2>  sin  B_ '  -  11 

U2  C  <  W2/U2)  sin  &2'  ~  13  +  )iU22  C-2Wsinfl2' +1*  W2  21 

U2  U2 

=  <W-?/Uo)  sin  Bo'  -  1 
-)4  ♦  >4  <W22/U22) 


'c  =  -  2  C  <Wo/Uo  sin  132'  -  1  3  A(14> 

<  1  -  W2/U2> 

Equations  A<10)  and  A<14)  are  valid  for  a  discreet  jet, 
or  for  distributed  gas  at  the  rim.  For  the  distributed  gas, 
it  is  possible  to  calculate  the  referred  power  as  a  function 
of  wheel  radius,  giving  specific  power,  P3.  Using  equation 
A<  10)  , 

Pa  =  P/mU2  =  C  <  W2/U2 )  sin  132'-  13  A<15) 

Therefore  equation  A(14>  becomes, 

"c  =  2  Pc. 

(W2/U2)2  ~  1 

For  sin  f32'  =  1*  implying  that  d2' s90  degrees, 

"c  3  2 

( W2/U2  >  *  1 
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therefore 


V©  =  W  sin  B  *  U 


A  ( 7 ) 


Substituting  equations  A(7)  into  ACS)  we  obtain, 
W*  =  P  =  -  U2  l  W  sin  B  U  3 


'x  -  JEL 

m 


-  U  2  I  W  sin  B  *  1  3 
U2 

U22  C  W  sin  <-B>  -  1  3 
U2 


A  <  8  ) 


Defining  the  power  conversion  efficiency  tc  as  the 
fraction  of  the  available  power  actually  converted  to  shaft 
work,  thus 

P/m _  A  ( 9> 


•c  = 


P/m  ♦  v22/2 


and  defining  B2'=  -B2»  the  relative  flow  angle  at  the  rim 

which  is  positive  in  the  direction  of  rotation,  then 

P  =  U2  [  W2  sin  B2'  -  1  3 .  A<10) 

m  U2 

From  the  geometry  shown  in  Figure  Al, 


so  that. 


V22  =  V©22*  (  W22  -£V02-U2]2) 


JiV22*  CU2  V ©2  -  HU22*  W223 


Using  equations  ACS)  and  AC10) 

-  U2  V ©2  =  -Wx  3  U2  I_W2  sin  fl2'  -  13 

U2 

Substituting  equation  AC12)  into  A(ll)  gives 


A(ll) 


A  C 12 ) 
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p 


w'  T 


A<2> 


The  work  per  unit  ease  flow  ia  given  by 

P  /  m  -  w'  Cri  V  9i  -  r2  V©2>  A<3) 

The  aubacript  1  denotea  entrance  station,  2  denotea  exit. 
Defining  the  rim  apeed  as  U  =  rw',  equation  A(3)  becomes 

P  /  m  =  Wx  =  Ui  V0i  -U2  V  02  A<4> 

A  general  velocity  diagram  for  a  rotor  ia  ahown  in  Figure 
Al. 

The  tangential  component  of  velocity  V0  ia  taken 
poaitive  in  the  direction  of  the  rim  velocity  U.  Figure  Al 
ahowa  all  poaitive  velocity  components  in  conventional 
cylindrical  coordinates. 

The  gas  enters  the  wheel  near  the  center,  effectively 
with  zero  moment  of  momentum  (about  the  axial  direction). 
Therefore  V 01  ■  O,  and  from  equation  A (4), 

«x  =  _P_  ■  -02  v®2  A (5) 

a 

For  Wx  to  be  poaitive,  it  ia  seen  that  V02  is  negative.  The 

exit  velocity  diagram  representing  the  flow  at  the  rim  of 

the  wheel  is  given  in  Figure  A2. 

From  the  geometry  ahown  in  Figure  Al,  in  general 

ain  B  =  V0  -  U  A(6> 

W 
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APPENDIX  A 


TORQUE  PRODUCED  BY  DISTRIBUTED  GAS  REACTION 

A.l  METHOD 

The  RDWE  concept  involves  an  unsteady  process.  The 
potential  performance  of  the  RDWE  can  be  understood, 
approximately,  by  examining  the  range  of  conditions  through 
which  the  unsteady  process  passes.  Thus  the  steady  flow 
performance  of  the  proposed  arrangement  needs  to  be 
examined.  We  first  attempt  to  predict  analytically  the 
power  or  torque  impulse  that  can  be  generated  from  the 
engine  under  steady  flow  outlet  conditions. 

The  passages  in  the  RDWE  exit  nearly  tangentially 
(Figure  1).  Since  this  condition  is  difficult  to  achieve  we 
examine  the  effect  of  the  exit  angle  on  the  amount  of  power 
that  can  be  generated  and  the  efficiency  that  can  be 
produced . 

The  basic  equation  for  relating  the  fluid  flow 
conditions  to  the  energy  transferred  to  the  shaft  is  Euler's 
Turbine  equation.  The  torque  of  the  fluid  acting  on  the 
rotor  is  given 

T  =  m  <  ri  V0i  -  r2  V&2>  A<1) 

If  the  rotor  has  a  rotational  speed  w',  the  shaft  power 
transferred  or  rate  of  work  transferred  from  the  fluid  to 
the  rotor  is  given  by. 


Figure  10.  Computational  Grid  for  Outflow 

Analysis 
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R0  ”4  Ri  *  3.2  No.  Bladaa  ■  60 

9.1.72521902  0  ■  80.82738092 

Condition  of  Tangancy  la 
sin  0  *  coa 

ain  0  *  0.159407755  coa  -  0.159409430 
chord  *  0.983540343  chord  taat  ■  0.983540073 
alpha  -  18.34514098 
radial  eolidity  ■  1.909859317 
thicknaaa  to  chord  ratio  ■  0.194933797 
radiua  of  caabar  circla  ■  0.758551242 
radiua  of  auction  circla  ■  0.575175292 
radiua  of  praaaura  circla  *  1.462802228 

Figura  7.  Final  Blading  Daaign. 


Figura  4.  CDTD  Flow  Genarator  -  Section  AA 


from  Figura  3. 


Figure  3.  Schematic  of  the  Proposed  High  Speed  CDTD. 


Thus  the  conversion  efficiency  is  unity  when  W2/U2  =  1 
but  this  is  a  condition  which,  from  equation  A(15)  gives 
zero  power  output. 

A. 2  RESULTS 

The  value  of  the  specific  power  and  conversion 
efficiency  (given  by  equation  AOS)  and  A(16))  as  a  function 
of  the  ratio  of  relative  flow  to  wheel  velocity,  for  various 
values  of  the  relative  exit  angle,  are  shown  in  Figure  A3 
and  Figure  A4  respectively. 

Figure  A3  shows  that  the  specific  power  increases 
linearly  as  the  exit  velocity  increases,  and  decreases  as 
the  relative  exit  angle  becomes  less  tangential.  The 
conversion  efficiency  is  seen  to  decrease  as  the  power 
increases  and  to  decrease  as  the  relative  flow  angle 
decreases  (less  tangential). 

The  results  shown  in  Figure  A3  and  Figure  A4  can  be 
used  to  gauge  the  performance  of  the  unsteady  RDWE.  In  the 
unsteady  expansion  process  the  value  of  the  relative  flow 
angle  (02>  will  be  constant,  being  determined  by  the  passage 
exit  geometry.  The  wheel  speed  will  be  approximately 
constant  and  only  the  value  of  the  relative  velocity  will 
change.  Thus  the  unsteady  process  at  different  times  will 
be  described  by  points  on  a  single  line  shown  on  Figure  A3 
and  A4,  depending  on  the  particular  value  of  f)2- 
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The  difficulty,  clearly,  ia  that  the  exit  velocity  must 
not  be  allowed  to  decrease  below  a  specific  factor  times  the 
wheel  speed  or  power  will  be  absorbed  rather  than  generated. 
Also  the  conversion  efficiency  will  vary  during  the 
expansion,  but  will  always  be  poor. 

This  implies  that  the  concept  of  the  RDWE  as  shown  in 
Figure  1  can  not  be  efficient.  The  device  would  require 
some  type  of  diffuser  that  would  convert  the  kinetic  energy 
leaving  the  rotor  into  static  pressure.  Unfortunately,  a 
variation  in  the  relative  velocity  magnitude  implies 
variation  in  both  magnitude  and  direction  of  the  exit 
absolute  velocity.  A  suitable  diffuser  for  such  a  flow  does 
not  exist.  Consequently,  other  geometries,  particularly 
those  which  would  allow  constant  flow  to  the  diffuser  by 
sequencing  the  detonation  wave  with  (partial)  ports  at  the 
rotor  exit,  should  examined. 


TABLE  A1 

Specific  Power  and  Power  Conversion  Efficienc 


Vs  velocity  Ratio  Data 


60  DEG 


70  DEG 


0 

-1 

2 

0 

-1 

2 

5 

-.56699 

1.51200 

.5 

-.53015 

1.4137 

0 

-.13397 

0) 

1.0 

- .06031 

a> 

5 

.29904 

.47846 

1.5 

. 40953 

.65526 

0 

. 73205 

.48803 

2.0 

.87938 

.58626 

5 

1.16506 

.44383 

2.5 

1.34923 

.51399 

5 

2.30109 

.36108 

3.5 

2.28893 

.40692 

5 

2.89712 

.30099 

4.5 

3.22862 

.33544 

5 

3.76314 

.25731 

5.5 

4.16831 

.28501 

5 

4.62917 

.22444 

6.5 

5.10800 

.24766 

5 

5.49519 

. 19892 

7.5 

6.04770 

.21892 

BO  DEG 


90  DEG 


0 

-1 

2 

0 

-  1 

2 

5 

-.50760 

1.35359 

.5 

-.5 

1.33333 

0 

- .01519 

CO 

1.0 

0 

1 

5 

.47721 

.76354 

1.5 

.5 

.80000 

0 

.96961 

.64641 

2.0 

1.0 

.66667 

5 

1.46202 

.55695 

2.5 

1.5 

.57143 

5 

2.44683 

.43499 

3.5 

2.5 

.44445 

5 

3.43164 

. 35653 

4.5 

3.5 

.36364 

5 

4.41644 

.30198 

5.5 

4.5 

.30769 

5 

5.40125 

.26188 

6.5 

5.5 

.26667 

5 

6.38606 

.23117 

7.5 

6.5 

.23529 

APPENDIX  B 


DOUBLE  CIRCULAR  ARC  BLADING  DESIGN 

B.l  METHOD 

A  computer  procedure  waa  developed  to  design  radial 
out-flow,  near  tangential  turning  vanes,  specifying  only  the 
outer  radius,  the  inner  radius,  and  the  number  of  blades. 

The  blade  spacing,  radial  solidity,  camber  angle,  suction 
and  pressure  surface  radii,  chord  length,  thickness  to  chord 
ratio,  and  efflux  angle  are  determined.  The  suction  surface 
is  designed  to  be  tangential  at  the  rim. 

A  Hewlett  Packard  9830A  computer  was  used.  The  blade 
shapes  were  generated  using  a  Hewlett  Packard  9862A  plotter. 

The  geometrical  approach  to  the  design  of  the  vanes  is 
illustrated  in  Figure  B1  and  much  of  the  notation  is  defined 
in  this  figure.  The  following  assumptions  were  made  as 
initial  constraints: 

1.  radius  at  exit  R0,  is  known 

2.  radius  at  entrance  Rj_,  is  known 

3.  number  of  blades  Z,  is  known 

4.  The  suction  side  is  tangent  to  the  outer 
radius  (point  C) . 

5.  The  camber  line  c',  is  tangential  at  the  inlet 
radius  (point  A) . 

6.  The  blade  subtends  an  angle  ^  at  the  center. 
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Initially,  a  procedure  was  developed  which  required 
specifying  the  caliber  angle  0,  R0,  Rj.,  in  addition  to  the 
nuaber  of  blades.  Since  the  choice  of  0  which  would  give 
acceptable  values  of  efflux  angle  and  blade  thickness  was 
not  obvious,  an  equation  for  the  condition  that  the  camber 
line  was  tangent  at  the  rim  (which  would  give  a  blade  of 
zero  thickness,  since  the  auction  surface  was  also  required 
to  be  tangent)  was  written.  The  equation  was  seen  to  have 
two  possible  roots.  The  selection  of  the  smaller 
(incorrect)  root  gave  positive  values  of  the  efflux  angle 
and  reasonable  values  of  blade  thickness.  The  program  which 
was  developed  provides  the  option  of  selecting  0,  or  of 
selecting  the  blade  which  results  from  finding  the  incorrect 
smaller  root  to  the  tangent  camber  line  condition. 

Described  here  are  the  equations  for  the  generation  of  the 
blade  shapes. 

Defining  peripheral  blade  spacing  s  as, 

s  =  2  n  R0  /  Z 

and  the  angle  subtended  by  the  blade  space  as 

fa  =  360  /  Z 

the  radial  chord  as 

f  ■  <Ro  -  Ri>  /  Ro 

then  the  radial  solidity  is  given  by 

(Tj-  =  (  Rq  -  Ri )  /  s  =  (O  / 
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Using  these  expressions  and  Figure  B1  the  following 


equations  can  be  written: 

GE  =  Rc  -  Rc  cos  (0/2) 

=  Rc  Cl  -  cos  (0/2) 

HE  =  c  r/2  *  Rc  £  1  -  cos  (0/23 
Where  T  is  the  thickness  to  chord  ratio.  Also 

HE  =  Ra  C  1  -  cos  ( 0s / 2  >  3 
To  solve  for  the  thickness  and  r  we  subtract  the 
equations  B(l)  and  B<2) 

HE  -  GE  =  t/2  *  t  c/2 

Using  Figure  Bl, 

FE  =  Rc  Cl  -  cos  <0/2)3  -  r  c/2 
=  Rp  Cl  -  cos  <  0p / 2 ) 3 

Subtracting  equation  B<5)  from  B(l), 

GE  -  FE  =  t/2  =  r  c/2 

The  chord  c,  is  given  by, 

c  *  2  Rp  sin  <0p/2> 

=  2  Ra  sin  <  0s/2) 

*  2  Rc  sin  <0  /2> 

All  equations  were  sade  dinensionless  by  dividing  by 
giving  for  exasple 

c/R0  =  c" 

c”  *  2  Rp.  sin  (0p/2> 

Ro 

*  2  Ra  sin  <0a/2) 

Ro 


B(l) 

B<2> 

B  <  3 ) 

B  <  4  ) 

B  ( 5 ) 
B(6) 

B  (  7  ) 

B  ( 8 ) 
B  ( 9 ) 
B  ( 10 ) 

t 

B  ( 1 1 ) 
B  <  12 ) 

B  ( 1 3 ) 
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B  <  14  ) 


=  2  Rc.  sin  <  0  /2> 

Ro 

For  the  geometry  we  seek,  given  R0,  and  the  camber 
angle  0  and  efflux  angle  a,  the  angle  ACO  (or  A)  is  given  by 

A  =  90  -  (  a  *  0/2 >  B (15) 

From  the  triangle  ACO 

f  *  90  -  (.a  *  0/2>  ♦  180  -  0/2  =  180 

resulting  in, 

90  ♦  -  0  ♦  a  B  ( 16 ) 

This  equation  fully  defines  a.  The  values  of  *{*  and  0 
are  needed  to  complete  the  description  of  the  camber  line. 

To  solve  for  the  radius  of  each  blade  arc  the  following 
relations  were  used. 

First  with 


and 


Aa 


Re-  -  Rr-  cos  0/2  =  HE 
Ro  Ro 


B  ( 17  > 


Ap  *  Rp  -  Rp  cos  < 0p/2)  =  FE  B(18) 

Ro  Ro 

Using  equation  B(3),  B(6>  and  B(7) 

Rc  £  1  -  cos  <0/2)3  -  Rp  t  1-cos  <0p/2>]  =  r  c/2  B<19) 

Also, 

Aa  -  T  c/2  -  Rp  £1  -  cos  <0p/2>]  =  T  c/2  B<20) 
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which  gives 


Ap  =  Aa  -  r  c  B ( 21 ) 

Using  this  relationship  between  Aa  and  Ap  we  can  solve  for 
Ra  and  Rp.  From  equation  B<2),  and  B<3>, 

As  =  Rs  Cl  -  cos  ( 08/2 )  B ( 22 ) 

or 

cos  <  0s/2 )  *  1  -  Aa/Ra 

Squaring  and  adding  together  equations  B<22),  and  B(10) 
we  obtain 

cos2  <  0s/2)  ♦  sin2  <0s/2>  =  <1  -  Aa/Ra>2  *  <c/2Rs>2 
or 

1  *  1  -  2<Aa/Ra>  ♦  <Aa/Ra)2  ♦  c2/4Ra2 
which  becomes 

Rs  =  HAs  ♦  c  2_1 _  B (23) 

8Aa 

or 

R«  *  J*Aa  ♦  c  2  1 B  ( 24 ) 

Ro  R©  QAa 

Using  equations  B(5),  B<6),  and  B(16)  a  similar  manipulation 
will  yield, 

Rp  »  U  Ap  ♦  c2  1  B ( 25 ) 

Ro  8Ap 

The  angles  0s  and  0p  are  found  using  equations  B(13), 
and  B(14)  since 

sin  <0s/2)  =  c/2Ra 
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or 


0s  *  2  arc  sin  <  c/Ro) 

2R3/R0 

and  similarly 

sin  <0p/2>  =  c/2Rp 

or 


B  ( 26) 


0p  *  2  arc  sin  <  c/R.-.)  B<27) 

2Rp/R0 

Ths  desired  condition  is  that  the  suction  surface  be 
tangent  to  the  outer  rim.  From  Figure  B1  this  is  given  by 
the  relation 


0s/2  =  a  +  0/2  B (28) 

T  is  not  an  Independent  parameter  since  it  can  be 

calculated  using  equation  BOS),  thus 

r  *  t  <Rjt  C  1  -  cos  (  0s/2)  3  -  Rr-  C  1  -  cos  (0/2)]  )  B(29) 

<c/2) 

To  insure  near  tangential  flow  at  the  exit  we  first 
seek  the  condition  of  a  tangent  camber  line  at  the  exit. 

This  implies  since  we  also  seek  a  tangent  suction  surface, 
that  the  blade  would  have  zero  thickness.  Using  the 
notation  in  Figure  B2  the  following  equations  apply: 

sin  (0/2)  =  L  /  c  BOO) 

and 

sin  =  L/  R0  B (31 ) 

Substituting  equation  BOO)  into  B(31)  gives. 


c  *  R«-%  sin  'V' 
sin  (0/2) 


B  ( 32 ) 
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and  dividing  equation  BOS)  by  B(39), 

Rn  ain  ^  =  tan  < 0/2)  B(39) 

Ro  coa  f  -  Ri 

This  equation  eatabliahea  the  relationship  between  ^ 
and  0.  Using  Figure  B3  the  relationship  of  c,  7^,  0  is 

established, 

cos  (0/2)  =  c/2  B ( 40 ) 

R  tan 

or 

c  *  2  Ro  tan  ^  cos  (0/2)  BC41) 

This  is  the  condition  for  the  camber  line  to  be  tangent  at 
the  exit. 

Using  equations  B(32)  and  B(41), 


65 


sin  =  2  tan  cos  (0/2)  B(42) 

sin ( 0/2) 


giving , 

cos  y  =  2  sin  (0/2)  cos  <0/2> 

~  2  i  H  sin  <  0/2  ♦  0/2)  sin  <0/2  -  0/2)  ] 

=  sin  0  B  <  43 ) 

Therefore  for  the  camber  line  tangency  condition 

cos  ^  =  sin  0  B ( 44 ) 


Using  equation  B<39>  and  B(44) 

2  sin  <  0/2)  cos  < 0/2) 

2  cos  2  <  0/2 ) 


we  obtain  the  following 
=  tan  ( 0/2) 


and 


2  sin  (0/2)  cos  (0/2) 
2  cos  2  ( 0/2) 


Ro  sin 
Ro  cos  y- 


V' 


Ri 


Using  the  general  tr ignometr leal  relationship, 

sin  2x  =  2  sin  x  cos  x 


we  can  manipulate  the  preceeding  equation  to  obtain 

3in  0  C  Ro  cos  Jt  -  R i  3  =  2  cos  2  (0/2)  B<45) 

R0  sin  'J* 

Using  equation  B(43)  this  becomes 

cos  y  C  Ro  cos  V  -  Rj_  ]  =2  cos2  (0/2)  B(46) 

Using  the  relationship 

cos2  x  =  sin2  x  *  cos  2x 

2  cos2  (0/2)  =  2  Csin2(0/2)  +  cos  0] 

=  2  I  1  *  cos  0  ] 

2 
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=  1  ♦  <1-  sin2 

Using  equation  B(43)  this  becomes, 

2  cos2  <  0/2)  =  1  +  <1  -  COS2  0)** 

=  1  +  sin  X*  B ( 47 ) 

This  gives 

R0  cos^  -  Ri  cos  y-  -  R0  sin  J-  *  RQ  sin2^  B(48) 

Writing  a-ain')^'  equation  (48)  becomes, 

R0  < 1  -  s2)  -  R i  <  l-s2>  *  =  R0  s  +  R0  s2  B(49> 

or 

Rja  [  <l-s2)  -  a2  -  a3  =  -<l-a2)** 

Ri 

( Rn ) 2  [  1  -  2s  -  3s2  ♦  4s3  ♦  4s*  ]  =  l-s2 

<Ri> 

Thus  equation  B<48)  gives  the  condition 

F  < s )  =  <Rj_)2  (l-a2)  -  1  ♦  s<2  ♦  3s  -  4s2  -  4s3  )  =  0  BOO) 
Ro 

The  condition  of  tangency  of  the  camber  line  is  given 
by  the  solution  of  equation  B<50)  for  s  between  0  and  1, 
since  s  =  sin  /Y  .  Using  Newton's  method  of  interation  the 
root  can  be  found  for  a  given  radius  ratio.  Figure  B4  shows 
F(s)  in  equation  BOO)  for  various  radius  ratios.  Figure  B4 
clearly  identifies  two  possible  roots.  The  larger  root  is 
the  correct  root,  giving  zero  blade  thickness.  The  smaller 
root  if  incorrect,  but  gives  positive  values  of  a  and  finite 
blade  thickness.  Emphasis  was  placed  on  finding  the  first 
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1 


Figure  B2 .  Geometrical  Deecription  for  Camber  Line 

Tangency  Condition. 
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TABLE 

Bl.  List  of  Computer  and  Text  Symbols  for 
Double  Circular  Arc  Blades  (con't) 

\CIRC" 

TEXT 

IDENTIFICATION 

T1 

r 

thickness  to  chord  ratio 

T2 

0S 

suction  angle 

T3 

0p 

pressure  angle 

T4 

0 

loop  parameter  to  draw  arcs 

T9 

- 

tolerance  on  Newtons  method 
of  integration 

XI 

Xds 

suction  center  x  coordinate 

X2 

XdP 

pressure  center  x  coordinate 

X3 

Xd 

camber  center  x  coordinate 

X4 

- 

variable  for  x  values  in 
continuous  loop 

X9 

- 

variable  for  x  values  in 
rotational  loop 

Y1 

Yds 

suction  center  y  coordinate 

Y2 

YdP 

pressure  center  y  coordinate 

Y3 

Yd 

camber  center  y  coordinate 

Y4 

- 

variable  for  y  values  in 
continuous  loop 

Y9 

- 

variables  for  y  values  in 
rotational  loop 

20 

Z 

number  of  blades 
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TABLE  B1 .  List  of  Computer  and  Text  Symbols  for 
Double  Circular  Arc  Blades  (con't) 


ACIRC" 

TEXT 

IDENTIFICATION 

02 

cos 

cos  ^ 

03 

Rc 

output  print  of  Rc 

04 

Ra 

output  print  of  Rs 

05 

sin 

for  tangent  test 

06 

COS  0 

for  tangent  test 

PI 

n 

n 

P4 

- 

constant,  cRo/2Rp 

R0 

Ro 

outer  radius 

R1 

Ri 

inner  radius 

R2 

density,  (R©  -  Ri>/R0 

R3 

Rc/Ro 

camber  radius,  dimensionless 

R4 

Ri/Ro 

inner  radius,  dimensionless 

R5 

R»/Ro 

suction  radius, dimensionless 

R6 

Rp/R0 

pressure  radius, dimensionless 

SO 

y. 

periodicity,  2nRo/Z 

SI 

<rr 

radial  solidity 

S2 

s 

spacing 

S3 

'f 

blade  curvature  with  respect 
to  camber  line 

S5 

y. 

periodicity,  360/Z 

TO 

0 

camber  angle 

TABLE  B1 .  List  of  Computer  and  Text  Symbols 
for  Double  Circular  Arc  Blades 


“ACIRC"  TEXT 


IDENTIFICATION 


AO  a 

A1  As 

A2  Ap 

CO  c“ 


Cl 

C2 

11 

12 
13 


c 

c 


K1 

K2 

K9 


L 

N  N 

NO 


N9 


01 


efflux  angle 

suction  line  to  chord  length 
pressure  line  to  chord  length 
c/R0  dimensionless  parameter 
for  chord  length 
chord  length 

chord  length  for  tangent  test 

Interval  ( 0/2  -  0s/2> 

Interval  ( 0/2  *  0s/2) 

Interval  (0/2  *  0s/2)  -  (0/2 

-  0s/2) 

constant  value 

Newtons  method,  y/y' 

initial  guess  for  Newton's 
method 

counter 

Number  of  blades  to  draw 

Number  of  blades  to  draw  for 
loop 

Number  of  points  for  interval 
of  arc 

Equals  1  for  solution  of  ,2 
for  input  y 
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910 

920 

Y  3  =  R  4 
REM** 

**DRA W  BLADE  S IJ R FACES**** 

930 

DISP 

"INPUT  THE  NUMBER  OF  BLADES 

TO  DRAW"; 

940 

WAIT 

1000 

950 

DISP 

..  « 

960 

DISP 

"STARTING  AT  H=l"» 

970 

INPUT 

N 

930 

REM— 

- EDIT  NEXT  LINE  FOR  PLOT 

RESOLUTION 

990 

N9  =  25 

1000  FOR  N 0=1  TO  N 
1010  1 1  =  <  T 0/2  > -  (  T 2/2 ) 

1020  12  =  <  T  0  /  2  >  +  <  T  2  2  > 

1030  1 3= < 1 2-1 1 >/N9 

1040  PEM***FOR  THE  SUCTION  SIDE*** 

1O50  FOR  T 4  = 1 1  TO  12  STEP  13 
106 0  X 4=X1-(R5* C 0 S ( T4  >  > 

1 0  7  0  V  4 = V 1  +  <  R  5  *  S I N  <  T  4  >  > 

1 030  X9=N4*C0S < < N0- 1 >  *S5)  +  < Y4*S I N< <  NO- 1 >  *S5 ) > 

1 0  9  0  Y 9 = Y 4  *  C  0 S  < < N  0-1 > *  S  5  > - < X 4  *  S I N < <  N  0 - 1 >  *  S  5 )  > 

1100  PLOT  X9 >  Y9 
1110  NEXT  T 4 
1120  PEN 

1130  F:EM****FOP  CAMBER  LINE**** 

1140  FOR  T4  =  0  TO  TO  STEP  TO  tl? 

1150  N4  =  N 3 - < R  3 *  C 0 S ■ T  4 > > 

1160  V 4  =  V 3  +  ( R 3 * 3 1 N  (  T 4  1 ) 

1170  X9=X4*C0S < < NO- 1 > *35 > + ( Y4*S I N< ■ NO- 1 > *35 > > 

1  1 30  Y9=Y4*C0S <  (NO- 1  >*S5  >- ( X4*S  I N <  ( NO- 1 ')  *35 )  > 

1190  PLOT  X9»Y9 
12O0  NEXT  T 4 
1210  PEN 

1220  FEM****FOR  PRESSURE  SIDE**** 

1230  1 4  =  T 0 / 2 - T 2 2 

1240  I5=T0/2+T3/2 

1250  1 6= U 5- 1 4, >/N9 

1260  FOR  T4= 1 4  TO  15  STEP  16 

1 270  X4  =  X2-< R6*COS ( T 4 ) ) 

1230  Y4  =  Y2  +  (P6*SItN:T4)  ) 

1290  : : 9  = : : 4 - C 0 3  < 1  N o -  1 >  * S 5 )  +  Y 4 * S I N  < < NO - 1 > *  S 5  > 

1  300  Y9  =  'i'4*C0S  <  ( NO-  1  >  *35 > -X4*S  I N (  ( NO-  1  ;•* 35  1 

1  305  PLOT  X9 >  Y9 

1310  NEXT  T 4 

1320  PEN 

1330  NEXT  NO 

13  35  FEM****■t*•‘'*■**-*•k***********■**  +  ■fc<■-■,■*-''■-■,■- 
1  340  END 

1  350  PEM---*-NENT 0N3  HILL  BE  USED  TO  CALCULATE*** 

1360  PEM*****ThNGENCY  ROOT  BETWEEN  0  S<  90  DECREES***-* 

1  370  FE M ***** S 0 L V INC  FOR  THE  1ST  ROOT**** - - 

1330  T9=0. O00O1 
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•420 

420 

440 

450 

460 

470 

4S0 

435 

490 

500 

510 

520 

530 

540 

550 

560 

570 

530 

590 

600 

610 

620 

630 

640 

650 

660 

670 

630 

690 

695 

700 

710 

720 

730 

740 

750 

760 

770 

730 

790 

800 

310 

820 

830 

840 

845 

350 

860 

870 

830 

890 

900 


C2=2*TAHS3*C0S  <  TQ/2  >  *R0 

REM******CQNCLUDES  CALCULATION  FOR  PSI  j  PHI > AND  C 

A  0  =  9  0  +  S 3-TO 

C0=C1-"R0 

R3=C0/(2*SIN(TO/2>) 

R4=R1/R0 
D0=ATH  (  R3/R4) 

REM****FOR  SUCTION  SURFACE  TANGENT  TO  OUTER  RIM* 
T2=<2*A0>-»  i0 
R  5  =  C  O  <  2  *  3 1 N  <  T  2  /  2  >  > 

T  1  =  (  ( R 5 * <  i  - C 0 S < T 2 2 >  >  ) - < R 3 * U - C 0 S<TQ / 2 >  >  >  > < C O / 2 
A 1  =  R  3  *  ( 1  -  C  0  S  <  T  0  2  >  >  +  C  O  *  T  1  2 
A2=A1-(C0*T 1 > 

IF  A 2 > O  THEN  640 
IF  L=1  THEN  620 

Ii ISP  "NO  SOLUTION-RI  IS  TOO  SMALL"! 

b i sp  "  " ; 

msp  ■ 

INPUT 
L=  1 

R 1  =  R 1  +  C  H  / 1 0  0  ) 

GOTO  20O 

R 6  =  (  A 2 / 2 )  +  ( C 0 1 2/  (  8*A2  > ) 

P  4  =  C  U  /  (  2  *  R  6  ) 

T  3  =  2  *  A  T  N  (  P  4  /  S  G  R  (  1  -  P  4 1 2  > ) 

0  3  -  R  3  *  R  0 
U  4  =  R  5  *  R  0 
Ql=R6*R0 
REM: 

REM*******F'RIHT  STATEMENTS- 


'H 

H 


INCREASES  RI  I .  E .  =25‘:  ENTER  H" ! 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT  " 
PRINT  " 
PRINT  " 
PRINT  " 
PRINT  " 
REM**** 


P0= " R0 " R I = " R 1 " NO .  EL ABES= " ZO 
PS  I  =  "S3 
PHI  = "  TO 

CONDITION  FOR  TANGENCY  IS" 

SIN  PHI =COS  PSI" 

SIN  PHI  ="05"  COS  PSI  ="06 
CH0RB=" Cl "CHORD  TEST="C2 
ALF'HA="A0 

RADIAL  SOLIDITY*” , R9 
THICKNESS  TO  CHORD  RATIO*"* T1 


RADIUS 

RADIUS 

RADIUS 

**END 


OF  CAMBER  CIRCLE  ="03 


OF  SUCTION  CIRC 
OF  PRESSURE  GIF 
PRINT  STATEMENTS 


REM******************: 
REM*t**C ENTERS  CAN  NON 


LE  =' 
C  L  E  = ' 


04 

01 


BE  DETERMINED** 


Xl=R3-< <R3*COS<T0/2) >-<R5*C0S<T2/2) ) 
Y 1  =R4  +  ( ( R3*C0S < T  0/2 ) > - < R5*C0S < T2/2) ) 
Y. 2  =  R  3  +  (  ( R 6 * C 0 S (T3/2 >  )-< R3*C0$ (TO/ 2 > ) 
Y2*R4-< ( R6*C0S ( T3  2 > > - ( R3*C0S ( TO  '2 ) ) 
::3=R3 


)*COSCTO/ 

3 1 N  C  T  0  3  > 
> *C  OS ( T  0  2 ) 
> *3 IN (TO  2) 


74 


■  ro 


8.2.4  Program  Listing 


Liatad  hara  ia  tha  progra*  liating, "ACIRC“  for 
tha  deaign  of  double  circular  arc  bladaa 


10  REM*******RCIRC  1  JUNE  1984********** 

20  REM*******PR0GRAM  TO  CRLCULRTE  8<  DRAW  RPC  BLADING******* 
30  REM******F0R  RN  INITIRLLY  OUTFLOW  CASCADE****** 

40  REM****** U SI H G  NEWTONS  METHOD  TO  SOLVE  FOR  THE******** 

50  REM****0F  PSI********* 

60  BEG 

70  SCALE  -1  .<  1 »  -1  >  1 

SO  D I S P  "SET  PLOTTER  700  BY  7<  Y) -< CONT >  "  5 
90  STOP 

95  REM***** INPUT  BRTR*** ********************* 

100  DISP  "ENTER  THE  OUTER  RADIUS"; 

110  INPUT  RO 

120  DISP  "ENTER  THE  INNER  RADIUS"! 

130  INPUT  R 1 

140  R 4  =  R 1 / R O 

150  DISP  "ENTER  THE  NUMBER  OF  BLADES “! 

160  INPUT  ZO 


170  L  =  0 

ISO  PEM********ENB  INPUT  DATA**** 

1S5  PEM*****CALCUL AT  IONS ******* ************** 

190  S5  =  360--'Z0 
200  PI =3. 141592654 
210  S2=2*P 1 *RO/ZO 
215  P9=(R0-R1  >'-'S2 
220  S1=<P0-R1 )/S2 
230  S0=32/R0 
240  R2=S0*S 1 

250  DISP  "WANT  PSI  CALCULATED  Y=1.0=N"» 

260  WAIT  1000 

270  INPUT  01 

2S0  IF  01=1  THEN  340 

290  DISP  "INPUT  PSI"! 

300  INPUT  S3 
310  02=C03S3 

320  T0s ATM <02/SQR< 1-02*2)) 

330  GOTO  3S0 
340  GOSUB  1350 

350  PEM*****NOW  DETERMINE  THE  VALUES  FOR  PSI >  PHI >  AND  C***** 
360  S3  =  ATN < K9  'SOP < 1 -K9 t 2 > ) 

370  T0=2*ATN<P0*SINS3/(R0*C0SS3-R1 ) ) 

380  05=S I NS3 
390  0 6  =  C 0 S T 0 

400  C1=R0*SINS3/SIN>:T0''2) 

405  REM************************************** 

410  PEM******TEST  FOR  T ANGFNCY***** 


73 


is  a  convex  arc.  To  remedy  this,  the  program  asks  the  user 
to  increase  the  dimension  on  the  inner  radius  by  a 
percentage.  The  program  will  use  this  value  to  solve  for 
the  required  inner  radius  to  achieve  a  suction  surface  that 
is  tangent  to  the  outer  rim. 

The  fifth  part  prints  the  output.  The  sixth  part 
draws  the  blade  surfaces. 

B.2.2  Required  Inputs 

The  following  inputs  are  requested: 

Outer  radius 
Inner  radius 
Total  number  of  Blades 

Yes  or  No  response  for  solving  for  the  value  of  *f- 
Percentage  increase  in  the  inner  radius  -  may 
be  required,  when  the  value  of  ¥  is  specified 
Number  of  Blades  requested  to  be  drawn 

B.2.3  Listed  Outputs 

The  following  outputs  are  listed: 

Outer  radius 
Inner  radius 
Number  of  Blades 
Angle 
Angle  0 


Results  of  test  for  tangency 


illustrates  the  computer  generated  shapes  and  passage  for  a 
60  blade  design. 

B.2  COMPUTER  PROGRAM  “ACIRC" 

B.2.1  Description 

A  computer  program,  "ACIRC"  was  written  to 
construct  double  circular  arc  blading  using  the  above 
procedure.  The  program  will  calculate  and  draw  blading  for 
a  initially  radial-outflow  using  Newtons  Method  to  solve  for 
a  suitable  value  of  ,  or  for  an  arbitrary  input  value  of  'f  . 
In  either  case  the  method  will  yield  a  suction  surface  which 
is  tangent  to  the  outer  rim. 

The  program  is  user  friendly  and  is  annotated  for 
ease  of  understanding.  The  program  is  divided  into  six 
major  parts. 

The  first  part  is  the  input  data  section.  The 
second  part  calculates  the  periodicity,  radial  solidity, 
spacing,  and  asks  the  user  if  the  value  of  ^  is  to  be 
calculated  or  to  be  specified. 

The  third  part  calculates  the  first  root  for 
tangency  conditions  at  the  exit.  No  input  is  required  for 
this  section.  The  fourth  part  calculates  the  parameters 
needed  for  the  suction  surface  to  be  tangent  to  the  outer 
rim.  Usually,  if  the  user  specified  a  value  of  ¥  there  is 
the  possiblity  of  the  distance  from  the  chord  to  the 
pressure  circle  being  negative.  If  this  occurs  the  result 
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u 

B 


r* 


X(j  =  Rq 

=  Rc^o 
yd  =  Ri 


=  Ri/Ro 

To  draw  successive  blades  on  a  rotor  the  program  must 

rotate  the  axes  about  the  rotor  center  after  each  blade  is 

drawn.  The  following  equations  were  used  (see  Figure  B5) . 

X  =  x  ♦  Y  sin  0 
cos  0  cos  0 

Y  ♦  x  sin  0  =  y 
cos  0  cos  0 


therefore, 

X  =  x  cos  0  *  y  sin  0 
Y  =  y  cos  0  -  x  sin  0 

Given  the  number  of  blades,  the  following  operations 
(for  all  three  arcs)  are  performed,  the  coordinates  are 
calculated  as 


x  =  *d,ds»dP  "  R  cos  0 
y  *  ydrdsrds  ♦  R  sin  0 

The  arcs  are  then  drawn  from  0  =  (  ( 0/2 >  -  ( 0s, p/2)  )  to 

((0/2)  «■  (0s, p/2)) 


using 


X  =  x  cos  (N)  360  ♦  y  sin  (N)  360 
Z  Z 


and 

Y  =  y  cos  (N)  360  -  x  sin  (N)  360 
Z  Z 


where  N  is  the  number  of  blades.  Figure  B6  and  B7 


t 


' 

i 


root.  The  value  found  will  give  the  selected  value  of  sin?1. 
Once  Y  is  found,  using  equation  B(39>  0  can  be  calculated 

and  c  is  obtained  using  equation  B(32). 

Using  equation  B<28)  and  B(9)  0s  and  Rs  can  be 
calculated  and  the  thickness  is  then  given  using  equation 
B ( 29 ) . 

A  test  is  made  to  assure  that  the  length  Ap  is 
positive.  If  Ap  is  negative  this  implies  a  convex  rather 
than  concave  pressure  side.  This  would  also  imply  that  the 
inner  radius  is  too  small  and  should  be  increased  to  produce 
concave  curvature  on  the  pressure  side.  Ap  is  found,  from 
equations  B<19)  and  B<16). 

The  alternate  procedure  did  not  require  the  camber  line 
to  be  tangent  at  the  exit.  In  this  case  the  value  of  Y'  is 
not  solved  for  but  is  left  as  an  input  parameter.  The 
values  of  0  and  c  are  then  found  directly  using  equation 
B ( 39  > . 

The  centers  of  curvature  are  given  by. 


*ds  =_R c.  ~ 

(  Rr 

cos 

(  0/2) 

~  Ra  COS 

<  0s/2) 

) 

cos 

<  0/2) 

. 

Ro 

Ro 

Ro 

• 

yds  =_R_l  * 

(  Rr- 

cos 

<  0/2) 

-  Rg  COS 

<  0s/2) 

) 

sin 

<  0/2) 

; 

Ro 

Ro 

Ro 

‘ 

XdP  =  Re.  * 

<  Rp 

cos 

( 0p/2) 

-  Be.  cos 

<  0/2) 

) 

cos 

(  0/2) 

Ro 

R0 

Ro 

• 

ydp  3  R_l  - 
Ro 

(  Rp. 
Ro 

cos 

<  0p/2) 

-  Be.  cos 
Ro 

<  0/2) 

) 

sin 

(  0/2) 

with 
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Figure  B3.  Continued  Geometrical  Description  for  Camber 

Line  Tangency  Condition 
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R0  =  4  Ri  *  3.2  No.  Blades  =  60 

Y'*  9.1.72521902  0  *  80.82738092 

Condition  of  Tangency  la 
ain  0  s  coa 

ain  0  =  0.159407755  coa  »  0.159409430 
chord  *  0.983540343  chord  taat  -  0.983540073 
alpha  «  18.34514098 
radial  aolidity  »  1.909859317 
thickness  to  chord  ratio  *  0.194933797 
radius  of  camber  circle  *  0.758551242 
radius  of  auction  circle  ■  0.575175292 
radius  of  pressure  circle  •  1.462802226 


Figure  B6.  Computer  Generated  Shape* 
■fox  s  60  D©ti9^ • 
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APPENDIX  C 


WEDGE-ARC  BLADING  DESIGN 

C.l  METHOD 

Because  of  the  limitations  seen  in  the  outlet  passage 
shapes  which  could  be  generated  using  purely  double  circular 
arc  blades,  symmetrical  linear  surface  segments  were  added 
in  an  attempt  to  achieve  a  convergent  nozzle  at  a  shallow 
angle  at  the  rim  exit. 

From  the  design  using  only  double  circular  arcs  it  was 
found  that  by  using  the  original  camber  line  and  the  suction 
line  to  describe  the  blade,  improved  convergency  at  the  exit 
was  obtained  (see  figure  Cl).  Removing  the  original 
pressure  line  would  cause  the  flow  to  be  at  a  small 
incidence  angle  at  the  inlet,  but  this  was  thought  to  be 
acceptable . 

By  allowing  the  flow  to  exit  over  a  straight  suction 
side  from  the  exit  of  the  passage,  better  control  of  the 
outlet  flow  angle  was  assured.  The  wedge  arc  geometry  is 
shown  in  figure  C2. 

Referring  to  figure  C2,  the  wedge  arc  blade  shape  was 
developed  to  have  the  following  features,  or  constraints: 

1.  The  pressure  side  is  a  circular  arc. 

2.  The  suction  side  between  point  H  and  I  is  a  circular 
arc,  between  straight  and  symmetrical  line  segments. 


3.  A  wedge  angle  (O  is  specified  at  the  leading 
and  trailing  edges. 

4.  Equal  leading  edge  and  trailing  edge  thicknesse  ( S) 
are  specified. 

5.  The  radius  to  the  suction  surface  at  the  point  of 
tangency  to  the  straight  segment,  (point  I) 
intersects  the  pressure  surface  of  the  adjacent  blade 
at  the  trailing  edge  (point  C'>,  or  downstream  of 
point  C',  depending  on  the  value  of  the  parameter 

(  C'  ) 

Referring  to  figure  C2  and  C3  the  following 
relationships  are  obtained: 


xp  - 

xc  “ 

S'  cos 

0 

C(l> 

yp  * 

yc  * 

S'  sin 

0 

C  (  2  ) 

where  5'  is  the  leading  edge  and  trailing  edge  thickness. 
The  value  of  S'  has  been  made  dimensionless  by  dividing  by 
Ro-  All  other  quantities  have  been  made  dimensionless  in 


the  procedures  given  in  appendix  B. 

xt  =  xc  -  S/2  cos  0  C (3) 

yt  *  yc  *  s/2  sin  0  C  ( 4 ) 

xe  *  -S  C(5> 

ye  ■  Ri  c ( 6 ) 

XL  =  -S/2  C(7> 

yL  =  Ri  C ( 8 ) 


88 


The  wedge  angle,  (  is  specified.  With  (  known,  the 
length  of  the  straight  section  <  X  )  from  the  exit  can  be 


found.  For  a  symmetric  blade,  the  length  <.  X  >  is  the  same 
at  the  inlet  and  the  exit.  Referring  to  figure  C2,  C3  and 
C4 

Jl  *  i  '<1  -  «')  C<9) 

where  &  -  Jl  ‘  when  the  point  I  lies  on  the  radius  C'D' 

(figure  C4> .  The  slope  of  the  line  C'F  (  p)  is  given  by, 

p  =  arc  tan  (  VP  -  Vr ' )  C(10) 

<  xf  -  xc'  > 

The  coordinates  of  C'  are  found  by  solving  for  the 
previous  blade  and  rotating  the  coordinates.  These  values 
are  placed  in  an  array  until  needed.  A  straight  line  is 
drawn  from  the  current  blade  to  the  previous  blade  as  is 
defined  as 

C'F  *  <  CyF  -  yc]2  ♦  CXF  ~  xc']2  >  *  C(ll> 

Referring  to  figure  C3,  the  slope  of  FG  is  given  by, 

A  =  to-  i'  -  i  1  C  ( 12 ) 

A  relationship  of  X  '»  C'F  and  the  line  FG  is  given  by, 

X'  -  C'F  cos  C' FG 

*  C'F  cos  C  a  -  ^  -  <  ♦  p]  C  <  13 ) 

From  equation  C<13>  the  coordinates  of  I  and  H  are  found, 

xi  =  xp  ~  i  cos  Cot  -  y*  -  (3  C ( 14 ) 
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K  "  A 


yi  =  yF  -  Q  ta  -  ^ 

-  n 

C(  15) 

xh  =  ~  I  S  ■»  JL  sin 

(3 

C  ( 16  ) 

yH  3  Ri/R0  ♦  &  cos 

( 

C  <  17 ) 

The 

centers  of  the  suction  side 

are  given  by 

either. 

xd'  3  xh  +  Rz  cos 

( 

C  ( 18 ) 

YD'  3  YH  ♦  Rz  sin 

( 

C  <  19) 

or. 

xd'  3  xi  ♦  R2  sin  I  a  - 

y  -  (3 

C  ( 20 ) 

yD'  3  yi  -  Rz  cos  [  0!  - 

y  -  (3 

C  ( 21 ) 

Fro> 

equations  C(17>,  C<19)  and 

C<11)  Rz  can 

be  found 

since 

xi  -  xh  =  Rz  (cos  (  - 

sin  ^  ) 

yielding. 

Rz  3  XT  -  XH 

cos  (  -  sin 

& 

C  ( 22 ) 

To  calculate  the  position  of  I  for  a  given  throat 
dlaenslon  <a')  at  the  exit,  further  manipulation  is 
required.  Referring  to  figure  C4,  where  “a"  denotes  the 
center  of  the  exit  throat, 

xa  =  xc'  +  a'  sin  a  C ( 23 ) 

and 

Ya  3  Yc'  ~  cos  a  C<24> 

Fro*  equations  C<21),  C<13)  and  C(15>, 

Rz  3  (  XF  -  i?  cos  a  S  ♦  i  sin  <  ) 

(cos  (  -  sin  ^) 

=  xf  ♦  S  ♦  Cain  (  -  cos  A ]  C(25) 

[cos  (-sin 
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Referring  to  figure  C4 

t  Rz  ♦  a']2  =  [  XD'_  xa  1 2  ♦  [  ya  -  YD'^2  C<26) 

Using  equation  C(20>  and  C(21)  the  above  equation  becomes, 

[  Rz  ♦  a']2  =  [xj-Rz  sin  A.  -  xa3  2  ♦  [ya-yj-Rz  cos  £1  2 
referring  to  equations  C<14),  and  C(15), 

t  Rz  ♦  a' l2  =  [xp  -  X  cos  A  ♦  Rz  sin  A  -xa32 

♦  £  ya  “  yF  +  lain  ^  -  Rz  cos  A  ]2 


which  reduces  to 


2Rza'  *  a '2  -  (xp  -  xa> 

2  ♦ 

<ya  - 

yf>2 

*  i2 

C  (  27  ) 

*2  1  cos  A 

<xa 

-  xF> 

♦  2  X 

sin  A 

<ya 

-  yF> 

♦  2RZ  cos  21 

<ya 

-  yp> 

♦  2RZ 

sin  21 

<  xp 

-  xa)  . 

Using  equation  C(25)  the  following  results  can  be 
obtained 


Ceos 


a'  2  * 


sin  A 


(xp-  xa)2  +  <ya-  yF>2  *  j.2  *  2  i  cos  A  (xa-XF> 


*2  X  Bin  A  <ya 

-yp) 

♦  2  f  XF  ♦  S  ♦ 

X  Csin 

€  - 

cos  A  1 

J  sin  A  (xp-xg) 

'  Ceos 

[  -  sin 

A  1 

j 

♦  2  /  XF  ♦  S  ♦ 

1  [sin 

C  - 

cos  A  ] 

)  cos  a  cya-ypi 

^  [cos 

(  -  sin 

A.  J 

4 

♦  2  /  XF  ♦  S  ♦ 

l[ain 

(  - 

cos  A  ]  1 

!  a' 

l  [  COB 

(  -  sin  A  ] 

J 

la' 2  *  C  coa 

(  -  sin 

» 

C<28> 

f  <*F  - 

xa>  2  ♦ 

<ya 

- 

yp>2 

•H*  N 

\»  2  1  cos  ^  ( xa 

-  xF 

CM 

♦ 

X  sinA  ( ya-yF)/ 

♦  2  xp  ♦  S  ♦  ACain  (  -  cos  A  ]_  sin  A(xf-x8) 

♦  2  xF  +  S  ♦  jf  tain  (  -  cos  A  1  coa  ^  (ya-yF) 

S  ♦  X Csin  <  - 


*2  xF  ♦ 


cos  A  3  a 


and  on  combining  terms, 

2 

Ceos  <  -  sin  ^  la'  =  C  cos  i  -  sin  A  3  C<29> 

♦  (2  Ceos  <  -  sin  ^  3  cos  A  (xa  -  xp) 

♦  2  Ceos  (  -  sin  A3  sin  A  <ya  -  yp)  )  Jl 
+  2  xp  sin  A(xp  ~  xa>  ♦  2  S  sinA(xp  -xa) 

♦  2  Ceos  (  -  sin  A  3  sin  A  <xp  -  xa) 

♦  2  XF  cos  A(ya  -yF>  ♦  2  S  co3A(ya  -  yp) 

♦  2  Ceos  (  -  sin  A  3  cosA(ya  -  yp) 

♦  2  xf  a'  ♦  2  S  a'  +  2  a'  Ceos  (  -  sin  A  3 

Ceos  (  -  sin  1  a'2  *  C  cos  (  -  sinA3  2  COO) 

♦  (  2  Ceos  (  -sin  A3  cos  A  <xa  -  xp) 

♦  2  Ceos  (  -  sin  A3  sin  A  <ya-  yp) 

♦  2  Ceos  (  -  sin  Al  sin  A<xf~  xa) 

♦  2  Ceos  (  -  sin  A  3  cos  A(ya  -yp) 

♦  2  a'  C  cos  {  -  sin  A3  } 

♦  2  xf  ainACxF  -  xa>  +  28  sin  A(xp  -  xa) 

♦  2  xp  cos  A<Ya  *  yF>  *  2  S  cos  A(ya  -  yp) 

♦  2  xp  a'  ♦  2  S  a' 

The  above  equation  is  quadratic  in  the  unknown  9  ,  and  a 
solution  can  be  written  explicitly.  The  value  of  Z  can  be 
found  using  equation  C(9).  The  result  is  as  follows. 

If 

B  =  C  2  Ceos  i  -  sin  A]  cosA(xa  -  xp)  COl) 
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♦  2 

Ceos 

(  -  sinA] 

sin  A ( ya 

-  y F> 

♦  2 

Ceos 

C  -  sinA] 

sin  A ( xp 

-  xa) 

♦  2 

Ceos 

(  -  sin  A ] 

cos  A (ya 

-  yF> 

♦  2  a  '  E  cos  (  -  sin  A  3  )  /  '^tcos  ( 


Let  C  = 


♦ 

then  our 


2  xp  sin  A  (xp  -xa>  +  25  sin  &(xp  -  xa) 

2  xp  cos  &<ya  -yF>  ■*•2  5  cos  A  (ya  ~  YF> 

2  xp  a'  *  2  6  a')/  '  2  Ceos  C  -  sin  ] 

quadratic  equation  is. 


0 


JJ'2Ccos  ( 
<1 


sin  _J  (l-(')2*  2B  ^/2Cc os 
(')  *  c  J^'2Ccos  (  -sin  ] 


-sin 


and  is  of  the  form, 

2 

0=  <!-(')  +  2  B  ( 1  -  (  '  )  ♦  C 


so  that 

(l-(')  =»  -2B/2  i  V(2B ? 2  -40 

2 

=  -  B  ±  V<  B  2-  C) 

The  following  possibilities  occur: 

1.  If  the  term  under  the  square  root  is  negative 
solution  does  not  exist. 

2.  If  the  term  under  the  square  root  is  positive 
the  value  of  <1  -  <')  is  between  0  and  1  then 


sinA  ] 

C  C  32 ) 


»  C ( 33 ) 

C  <  34 ) 

0(35) 


0(36) 

then  a 

and 


solution  exists',  the  solution  is 


3.  If  the  value  of  (1  -  (1>  obtained  in  step  (2)  is 

outside  the  range  of  0  to  1  then  the  solution  is, 

<  1  -  tl)  =  -B  -  )(B2-  C> 

4.  If  the  value  of  (1  -(1)  obtained  in  step  (3)  is 
outside  the  range  of  0  to  1  then  no  solution  exists. 


Once  ('  has  been  determined  and  all  other  values  have 
been  found,  the  design  is  complete  and  the  blades  can  be 
drawn . 

Referring  to  figure  C2,  to  draw  the  radius  on  the 
leading  edge  the  points 


X 

=  xl  ♦ 

S/2 

cos 

d 

C  <  37 ) 

y 

=  yL  - 

S/2 

sin 

6 

C  ( 38  ) 

are  plotted  over  the  interval  of  0  between  0  and  180 
degrees.  The  radius  on  the  trailing  edge  is  given  by. 


x  *  xt  - 

S/2 

cos 

0 

C  <  39 ) 

y  =  yT  ♦ 

S/2 

sin 

0 

C  ( 40 ) 

plotted  over 

the  interval  of  0 

between 

0  and 

0+180. 

To  draw 

successive  blades 

on 

the 

rotor 

the  procedure  is 

to  rotate  the  axes  after  each  blade  is  drawn.  The  procedure 
is  described  by  equation  B(55)  and  B(56). 
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C .  2  COMPUTER  PROGRAM  "AWEDG" 

C.2.1  Description 

The  wedge  arc  blade  design  program,  "AWEDG" 
calculates  and  draws  wedge  arc  blading  for  an  initially 
radial-outflow  rotor.  The  difference  between  "AWEDG"  and 
"ACIRC"  is  not  only  the  specified  input  parameters;  "AWEDG" 
calculates  the  straight  section  needed  on  the  suction  side 
to  constrict  the  flow  to  a  know  dimension  at  the  trailing 
edge  tip  of  the  adjacent  blade.  The  program,  for  simplifi¬ 
cation,  is  divided  into  seven  parts  and  each  part  will  be 
described . 

The  first  part  (lines  1  to  460)  contains  the 
values  of  all  specified  parameters.  On  running,  it  also 
allows  the  user  to  specify  the  angles  and  0  to  determine 

early  in  the  program  what  the  value  of  a  will  be.  If  the 
user  wishes  to  change  the  value  of  a  he  can  do  so.  Chord 
length,  camber  radius  and  a  are  calculated  in  this  section. 

The  second  part  (lines  470  to  640)  prints  some 
of  the  output  parameters  and  calculates  the  center 
coordinates  of  the  canter  line. 

In  order  to  calculate  the  outlet  passage  size 
the  preceding  passage  must  be  analyzed.  The  program  does 
this  in  part  three  (line  650  to  800)  by  calculating  the 
position  of  the  trailing  edge  tip  of  the  preceding  pressure 
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line 


Point  C'  on  figure  C3  is  calculated  by  rotation  of 


axes . 

In  part  four  (lines  810  to  1020)  the  position  of 
points  C,  F,  T,  E,  and  L,  shown  in  figure  C3,  are 
calculated. 

Part  five  (beginning  on  line  1030)  starts  the 
plotting  routine.  The  number  of  the  blades  requested  to  be 
drawn  are  specified.  The  pressure  side  radius  is  calculated 
and  is  drawn.  Calculations  for  the  straight  line  segments 
are  performed  (lines  1340  to  1570).  Based  on  the  determined 
value  of  the  straight  segment  the  coordinates  of  point  I  are 
dervied  (lines  1830  to  2090). 

In  part  six  the  radius  of  the  suction  side  is 
drawn.  Plotting  for  the  leading  and  trailing  edge  radius 
are  also  accomplished  in  part  six. 

Part  seven  (2100  to  2140)  concludes  the 
procedure  by  printing  the  value  of  the  throat  diameter, 
length  of  the  straight  line  segmnent  and  radii  of  pressure 
and  suction  circle. 

C.2.2  Required  Inputs 

The  following  inputs  are  required: 

Outer  radius 
Inner  radius 
Total  number  of  blades 
angle  ? 


angle  0  (must  be  consistent  with ^ -see 
discussion) 

Diameter  of  the  circle  at  the  passage  exit 
Leading  edge  and  Trailing  Edge  thickness 
Leading  edge  and  Trailing  Edge  wedge  angle 
Number  of  blades  requested  to  be  drawn 

C.2.3  Listed  Outputs 

The  following  quantities  are  output: 

Outer  radius 
Inner  radius 
Number  of  blades 
Angle  */- 
Angle  0 
Angle  a 

Radial  solidity 
Chord  length 

Thickness  to  chord  ratio 

Leading  edge  and  trailing  edge  thickness 

Leading  edge  and  trailing  edge  wedge  angle 

Throat  width  at  exit 

Radius  of  the  pressure  surface 

Radius  of  the  suction  surface 

Length  of  the  straight  section  on  the 

auction  side. 
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C.2.4 


Li at ad  hara  ia  tha  program  "AWEDG"  for  tha 
daaign  of  wadga  ahapad  ara  bladas.  Table  C2  idantlfiaa  tha 
ayabola  uaed  in  the  program  with  tha  defined  paramatera. 
Figure  C6  and  C7  illuatrataa  computer  generated  shapea  for 
an  outer  radiua  of  4  inches,  an  inner  radiua  of  3.2  inches 
and  60  blades.  The  value  of  ¥  and  0  are  specified 
parameters  to  limit  the  efflux  angle  a  to  between  8  and  9 
degrees. 


10  REM*******fiWEDG  1  AUGUST  1984********** 

20  REM***PROGRAM  TO  CALCULATE  &  DRAW  ARC  BLADING*** 

30  REH***ARCS  AND  STRAIGHT  BACKING  FOR  AN*** 

35  REM*** I  NIT  I  ALLY  OUTFLOW  CASCADE***** 

40  DEG 

50  DIM  GC  10Q] >  HE  1003*  UC  100  3 j  VC  100] 

60  SCALE  - 1 >  1  *  - 1  *  1 

70  DISP  "SET  PLOTTER  7  00  BY  7< Y)-<CONT> ”  5 
SO  STOP 

100  DISP  "ENTER  THE  OUTER  RADIUS" 5  * 

110  INPUT  R0 

120  DISP  "ENTER  THE  INNER  RADIUS"; 

130  INPUT  R 1 
140  R4*R1/R0 

150  DISP  "ENTER  THE  NUMBER  OF  BLADES" » 

160  INPUT  Z0 
170  L=0 
180  S5  =  360/'Z0 
190  P 1 =  3 . 141592654 
200  S2=2*F1*R0/Z0 
210  R9=(R0-R1 >/S2 
220  S1*<R0-R1)/S2 
2  J  0  _*  0  =  '*  2  /  R  0 
240  R2=S0*S1 

250  DISP  "ENTER  PSI-<  BTWN  LE  &  TE  0-360 /Z"! 

260  INPUT  S3 

270  DISP  "ENTER  PHI-PRESSURE  ANGLE- 'CAMBER< " » 

280  INPUT  T0 
290  A0=90+S3-T0 

300  PRINT  " ALHPA  IS  DESIRED  BWTN  8-9  DEGREES.  THE  VALUE" 
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p 


310  PRINT  "OF  PSI&PHI  CHOSEN  HILL  YIELD  AN 
320  PRINT  "ALPHA* "AO "DEGREES" 

330  PRINT  "PSI  =  " S3" AND  F'HI  =  "T0 

340  DISP  "NISH  TO  CHANGE  PSI&PHI  Y*1»N=0"! 

350  INPUT  01 

360  IF  01=1  THEN  250 

370  DISP  "ENTER  THE  LE  $<  TE  THICKNESS"  5 
3S0  INPUT  Q 1 


390  M0=Q1/R0 

4O0  DISP  "ENTER  LE  S<  TE  HEDGE  ANGLE "5 
410  INPUT  Ml 


420  Cl*  R  0  *  S I N  S  3  •••'  S I N  ( T  0  /  2 ; 

430  C0=C1/R0 
440  R  3  =  C 0 / < 2  *  S I N < T  0 / 2 ) > 

450  R4=R1/R0 
460  O3=R3*R0 
470  R  E  M  *  *  *  i 

4S0  REM*******F'R  I  NT  STATEMENTS 
490  PRINT  " R0= " R0 " R I = " R 1 " NO .  BLADES="Z0 
5O0  PRINT  "PSI="S3 
510  PRINT  "PHI ="T0 


520  PRINT  "CHORD* "Cl 

530  PRINT  " ALPHA* “AO 

540  PRINT  "RADIAL  SGLIDITY="R9 

550  PRINT  "LE  &  TE  THICKNESS* "01 

560  PRINT  "LE  &  TE  WEDGE  ANGLE* "Ml 

578  DISP  "ENTER  THE  THROAT  LENGTH  A* "5 

580.  INPUT  A5 

590  R|M******END  PRINT  STATEMENTS***** 

610  REM****CENTERS  CAN  NOW  BE  DETERMINED 
620  X3=R3 
630  Y3=R4 


640  REM**** DR AW  BLADE  SURFACES**** 

650  REM***FIRST  CALCULATE  COORDINATE  OF  C»  ON  BLADE  NO. 

660  N9=16 

670  NS =7 

6S0  N 1 =Z0 

690  1=1 


700  FOP  T 4=0  TO  TO  STEP  TO  N9 
710  X4*X3- < R 3 * C 0 S  < T4>  > 

720  Y4  =  Y3+  <  R3*S I N <  T4 ) > 

730  K9*X4*C03 < < N 1 - 1 > *S5  >  +  < Y4*S 
740  Y9=Y4*C0S < < N 1  - 1 > *S5 > - < X4*S 
750  PLOT  X9.Y9 
760  UC  I  ]=X9 


IN<<N1-1) 
I N  <  <  N 1  - 1  ) 


770  VC  I  ]=Y9 


7S0  1=1+1 
790  NEXT  T 4 
BOO  PEN 


>> 


Z-l* 
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305  REM***F  &  C  ON  PRESSURE  LINE*** 

310  REM***SECONO  CALCULATED  THE  COORDINATES  OF*** 
820  N0= 1 


8301-1 

840  FOR  T4=0  TO  TO  STEP  T0.-N9 
850  X 4 = X  3  - < R3*C0S ( T 4  > ) 

S60  Y4=Y3  +  <  R3*S I N < T4  >  > 

870  X 9 = X  ■ 4  *  C 0 S  < ( N 9 - 1 >  *  S 5  >  +  <  Y  4  *  S I N < 

8  8  0  Y  9 = Y 4  *  C 0  S  < < N  0-1 > *  S 5  > - ( X 4  *  S I N  < 

390  GE I 3=X9 

900  HE  I 3=Y9 

910  1=1+1 

920  NEXT  T 4 

930  M2=GC  N9+1 3 

949  M 3 = M 2 - M 9  *  C 0  S  T  9 

950  M4=M2-(M0''2>*COSTO 
960  L2=HC  N9  +  1 3 

970  L3=L2+M0*SINTO 
930  L4=L2+  ( N0-'"2)  *S  I NT0 
999  M5=-M9 


H0-1 >*S5> ) 
H  0  - 1 )  *  S  5  )  > 


1 0  9  0  M  6 = -  (  M  0  /  2  > 

1G10  L5=E1'"R0 
1920  L6=R1/R0 

1030  REM*****PLOTT I NG  SECTION***************** 

1040  DISP  "INPUT  THE  NUMBER  OF  BLADES  TO  DRAW " 5 
1950  INPUT  N 

I960  REM - EDIT  NEXT  LINE  FOR  PLOT  RESOLUTION - 

107O  N9-16 

1080  FOR  N0= 1  TO  N 

1090  REM****FOR  PRESSURE  LINE**** 

1190  1=1 

1110  FOR  T4=0  TO  TO  STEP  T0/N9 
1120  X4=X3- <  R3*C0S < T4 ) ) 

1139  Y4=Y3+  < R3*S I N < T4 ) > 

1140  X9=X4*C0S < i NO- 1 > *S5 >  +  < Y4*S I N < <N0-1)*S5) > 

1150  Y9=Y4*C0S ( <N0- 1 >*S5> - <X4*S1N< < NO- 1 >  *S5 > > 

1169  GE I 3=X9 

1170  HE  I 3=Y9 
1180  1=1+1 
1190  PLOT  X9«Y9 
1200  NEXT  T4 
1210  PEN 

1220  REM***SUCT I  ON  SIDE  LINE***** 

1230  REM***VALUE  OF  L  SOLVED  FOR  MILL  BE  GENERATED*** 
1240  REM***TO  GIVE  THE  DESIRED  THROAT  AT  EX  IT******* 

1259  E1=L3-VE N9+1 3 

1260  E2=M3-UE  N9+ 1 3 
1270  L0=RTN<-E1/E2> 

1230  C3=SQR<Elt2+E2t2) 


1290 

1300 

1310 

1315 

1320 

1330 

1340 

1350 

1360 

1370 

1 3S0 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1430 

1490 

15O0 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1530 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1630 

1690 

1700 

1710 

1720 

1730 

1740 

1750 


B 1 =R0-33-M 1 

B9=C3*C0S< A0-S3+L0-M 1 ) 

REM***CALCULATIGN  OF  BLADE  SUCTION  SURFACE*** 
REM***F0R  A  GIVEN  THROAT  SIZE  AT  EX  IT*** 

REM**#A4  IS  THE  RADIUS  AT  THE  THROAT  A5-'2*** 
A4=A5-'<2*R0> 

X6=UC  N9+ 1  ]+A4*SINA0 
Y6=VC N9+1  ]-A4*COSA0 

Z5=  < S I NB 1  *  < Y6-L3  >  +  COSB 1  * < X6-M3 )  )  *  <  C OSM 1 -S I NB 1 > 

Z  6  =  < S IND1* < M 3 - X  6 > +  0 0  S  Ii  1  *  < Y  6  -  L  3  >  -  A  4  >  *  < S I N  M 1 - C  0 S B 1 ) 
O 1  =  CZ5+Z6 > /  (  C3*  <  COSM 1  -S I  NB  1  >  ) 

Z7  =  < 2* ( M3+M0 ) * ( S I NB 1  *  < M3-X6 ) +COSD 1  * ( Y6-L3  > - A  4  > ) 
ZS=< < M3-X6  > t2  +  ( Y6-L3 ) T2-A4 12  >  * <  COSM 1  —  SI NB 1 > 
Q2»<Z7+Z3>/<C3t2#<C0SMl-SINBl ) ) 

03=0 1  t'2-Q2 

IF  03 <0  THEN  1560 

0  4  =  ~  0 1  +  S  U  R  (  0  3 ) 

IF  0<O4< 1  THEN  1470 
IF  04 <0  OR  04 >1  THEN  1500 
REM***FIRST  SOLUTION  FOR  E«<l-04)*** 

T9= 1 -04 
GOTO  1550 
Q5=-01-S0R(Q3) 

IF  0< Q5< 1  THEN  1530 
IF  Q5< 0  OR  05 >1  THEN  1560 
REM##*SECOND  SOLUTION  FOR  E=<l-Q5>*** 

T9= 1-05 


1560 


GOTO  1570 

PRINT  "NO  SOLUTION  EXISTS" 

B3=B9*  < 1 -T9  > 

M7=M3-B3*C0SD1 
L7=L3-B3*S I NB 1 
M3  =  - <  M0  +  B3*S I NM 1 ) 

L  3  =  R 1  /  R  0 + B  3  *■  C  0  S  M 1 

REM***RAB I  US  OF  SUCTION  SIBE*** 

R  7  =  < M 7 - M  3 > / < C 0 S  M 1 - S INB1 ) 

REM***CEMTERS  OF  SUCTION  RADIUS*** 

X  5 = M  3  +  R  7  *  C  0  S  M 1 
Y5=L3+R7*SINM1 

REM***F'LOT  ROUTINE  TO  DRAW  SUCTION  SIDE** 
REM***F'LOT  RADIUS  ON  LE  TIP*** 

FOR  T4=0  TO  130  STEP  130.NS 
X4=f16+  ( M0/2  >  *COST  4 
Y4=L6-  <  M0/2  *S  I  NT 4 

X9  =  X4*C0S ( ( N  O - 1 >  *  S  5 >  +  <  Y  4  *  S I N < <N0-1)*S5>) 
Y9=Y4*C0S < < NO- 1 ) *S5 > - C X4*S I N < < NO- 1 ) *S5  >  > 
PLOT  X9>  Y9 
NEXT  T4 
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1760 
1770 
1780 
1790 
1800 
1810 
1820 
1 830 
1840 
1850 
I860 
1870 
1380 
1390 
1900 
1910 
1920 
1 930 
1940 
1950 
I960 
1970 
1 980 
1990 
2000 
2010 
2020 
2O30 
2040 
2050 
2060 
2070 
2080 
2990 
2100 
2110 
2120 
2138 
2140 


X4-MS 

Y4=L3 

X 9 = X 4*C0 8 < < N 0 -1>* 8  5  >  +  C Y 4  * 
Y9=Y4*C0S C ( N0- 1 )  *85 > - < X4* 
PLOT  X9»Y9 

REM***PLQT  SUCTION  RADIUS 

1 1  = — M 1 

I2=2*M1+T0 

I3=< 12-11 )■••  N9 

FOR  T4= 1 1  TO  12  STEP  13 

X  4 = X  5  -  R  7  *  C  0  8  T  4 

Y4=Y5+R7*SINT4 

X9=X4*C0S  < C NO- 1 ) *85 )  +  <  Y4* 

Y9=Y4*C08<  CNO-l >*85)-<X4* 

PLOT  X9»Y9 

NEXT  T 4 

X4=M3 

Y4=L3 

X  9 = X  4  *  C  0  8  <  ( N  0  - 1  >  *  S  5  >  +  <  Y  4  •> 
Y  9  =  Y  4  *  C  0  8  <  <  N  O  - 1  )  *  8  5  )  -  <  X  4  * 
PLOT  X9..Y9 

REM***F'LOT  RADIUS  ON  THE 

1 1  =T0 

I2=TO+180 

I3=( 12-11 )/N3 

FOR  T4= 1 1  TO  12  STEP  13 

X4=M4-  (.  M/2  >  *  COST  4 


SIN< CNO-l )*S5) ) 
SI N<  <N0- 1 ) *S5 ) > 


8 1 N  C  < N0- 1  ) 
S I N  C  <:  NO- 1  > 


•85)) 
■85)  > 


•3 1 H  <  < NO  - 1 )  ■ 
•8  INC  (NO-l >  • 

TE  TIP*** 


Y4=L4+<M0/2)*SINT4 
X  9 = X  4  *  C  0  8  <  <  N  O  - 1  )  *  S  5  > 
Y  9 = Y  4  *  C  0  8  <  C  N  0  - 1 )  *  S  5  ) 
PLOT  X9>  Y9 
NEXT  T 4 
PEN 

NEXT  NO 
PRINT  "H*=" A5 
PRINT  "L="B8 
PRINT  "RADIUS  OF  PRE 
PRINT  "RADIUS  OF  SUC 
END 


+  <Y4: 
-<X4t 


8  INC  CN0-1 ) 
8 INC.  CN0-1 ) 


*85)) 

*85)) 


8 SURE  CIRCLE  =" 
TION  CIRCLE-" <R 


7*R0) 


TABLE  Cl.  List  of  Computer  and  Text  Symbols 
for  Wedge-Arc  Blade. 


AWEDG"  TEXT 


IDENTIFICATION 


AO  a 

A4 

AS  2*a' 

B8  J? 

B9  r 

CO  c" 


Cl  c 

D1 

El 

E2 

G  CIl 

H  [II 

I 

II 
12 


efflux  angle 

throat  dimension  (radius  of 
outlet  circle) 

throat  diameter  (diameter  of 
exit  circle) 

length  of  straight  section 

length  of  straight  section  to 
be  tangent  to  outlet  circle 

c/Ro  dimensionless, 
chord  length 

chord  length 

constant  angle 

constant  value 

constant  value 

array  for  current  pressure 
side  x  coordinates 

array  for  current  pressure 
side  y  coordinates 

counter 

Interval  variable 
Interval  variable 


13 


Interval  variable 


TABLE  Cl.  List  of  Computer  and  Text  Symbols  for 
Wedge-Arc  Blade  <con't>. 


AWEDGE" 

TEXT 

IDENTIFICATION 

L2 

yc 

final  value  on  present 
pressure  line,  y  coordinate 

L3 

yf 

final  value  on  present 
pressure  line  plus  blade 
thickness,  y  coordinate 

L4 

yt 

center  point  of  trailing 
edge,  y  coordinate 

L5 

yE 

Leading  edge  initial  suction 
side,  y  coordinate 

LG 

yi 

point  of  tangency  for 
suction  side,  y  coordinate 

L7 

yi 

point  of  tangency  for 
suction  side,  y  coordinate 

L8 

yH 

initial  point  for  suction 
side  curve,  y  coordinate 

HO 

s 

Leading  and  trailing  edge 
thickness 

Ml 

€ 

Leading  and  trailing  edge 
wedge  angle 

M2 

xc 

final  value  on  present 
pressure  line,  x  coordinate 

M3 

xf 

final  value  on  present 
pressure  line  plus  blade 
thickness,  x  coordinate 

M4 

*t 

center  point  of  trailing 
edge,  x  coordinate 

M5 

xe 

Leading  edge  initial  suction 
aide,  x  coordinate 

TABLE  Cl.  List  of  Computer  and  Text  Symbols  for 
Wedge-Arc  Blades  <con't>. 

“AWEDGE"  TEXT  IDENTIFICATION 

M6  xl  center  point  of  leading  edge 

M7  xj  point  of  tangency  for  suction 

side,  x  coordinate 

M8  xh  initial  point  for  suction 

aide  curve,  x  coordinate 


N 

N 

Number  of  blades  to  draw 

NO 

- 

Number  of  blades  to  draw  for 
loop 

N9 

- 

Number  of  points  for  interval 
of  arc 

01 

- 

Equals  1  to  change  ^  and  0, 

0  for  no 

03 

Rc 

output  print  of  Rc 

04 

Ra 

output  print  of  Rs 

PI 

IT 

n 

Q1 

- 

constant 

02 

- 

constant 

03 

— 

constant,  value  under  square 
root  term  for  solving  for 
<1 

Q4 

- 

same  as  03 

05 

- 

asms  as  03 

RO 

Ro 

outer  radius 

TABLE  Cl 


List  of  Computer  and  Text  Symbols  for 
Wedge-Arc  Blades  (con't) . 


AWEDGE" 

TEXT 

IDENTIFICATION 

R1 

Ri 

inner  radius 

R2 

e 

radial  chord,  <R0  -  Ri>/R0 

R3 

Rc^Ro 

camber  line  radius. 

dimensionless 

R4 

Ri/Ro 

inner  radius,  dimensionless 

R6 

Rp/R0 

pressure  radius, dimensionless 

SO 

s 

peripheral  spacing,  2nRo/Z 

SI 

(rr 

radial  solidity 

S2 

B/R0 

dimensionless  peripheral 

spacing 

S3 

f 

blade  curvature  with  respect 
to  center  of  rotor 

S5 

u 

spacing  angle,  360/Z 

TO 

0 

camber  angle 

T4 

- 

loop  parameter  to  draw  arcs 

T9 

- 

value  of  1- (' 

X2 

*dP 

pressure  surface  center,  x 

coordinate 

X3 

Xd 

camber  line  center,  x 

coordinate 

X4 

- 

variable  for  x  values  in 
continuous  loop 

X9 

- 

variable  for  x  values  in 
rotational  loop 
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TABLE  Cl.  List  of  Computer  and  Text  Symbols  for 
Wedge-Arc  Blades  (con't) . 


AWEDGE" 

TEXT 

IDENTIFICATION 

Y2 

ydP 

pressure  center  y  coordinate 

Y3 

yd 

camber  center  y  coordinate 

Y4 

- 

variable  for  y  values  in 
continuous  loop 

Y9 

- 

variables  for  y  values  in 
rotational  loop 

20 

2 

number  of  blades 

V 


Flgura  C2. 


W«dga-Arc  Shap«d 


Blading  G«o**try. 
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Not*:  C'D'i*  not  noriil  to  the  pr«*»ur«  surface 


Figure  C4.  Geometrical  Relationahip  Between  Throat 
Dieenaion  and  Suction  Surface  Radius 


R0  =  4  Ri  »  3.2  No.  of  Blades  *  60 
/  =  9.17  a  =  90.68 
chord  =  0.896198294 
alpha  =  8.49 

radial  solidity  =  1.909859317 
L.E.  &  T.E.  thickness  =  0.02 
L.E.  &  T.E.  wedge  angle  =  5 
A»  =  0.05225 
L  =  0.037306915 

radius  of  pressure  circle  =  0.629980625 
radius  of  suction  circle  =  0.476827358 


Figure  C5.  Computer  Generated  Shapes  for  a  Cascade  of  Wedge 
Arc  Shaped  Baldes  with  60  Blades. 
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